/ 



f (NASA-CE-151743) '^SHaTTLE KO-BAND SIGUAL 
DESIGN STUDY Final Report (LinCom Corp,, 

Pasadena, Calif.) 217 p HC AlO/HF A’01 
I CSCL 22B 

, G3/16 


N78-27162 ( 

\ 


Unci as 
25164 




PO Box2-793D, Pasadena, Calif 91105 



‘a£in(^c 


om 


FINAL REPORT 


SHUTTLE KU-BAND SIGNAL DESIGN STUDY 


PREPARED FOR 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
JOHNSON SPACE CENTER 
Houston, TX 77058 

Technical Monitor: William Teasdale 

Prepared Under Contract No. NAS 9-14846 


PREPARED BY 

William C. Lindsey 
Walter R. Braun 
Terri e M. McKenzie 




LINCOM CORPORATION 
P.O. Box 2793D 
Pasadena, CA 91105 




MAY 15, 1978 


TR-7805-0476 




jCinCc 


om 


ACKNOWLEDGEMENT 


The authors wish' to thank Dr. Bart Batson and William 
Teasdale of the Johnson Space Center for their support and for 
participating in various technical discussions during the course 
of this study. 


csClnG^i 


‘<jCinC^c 


om 


TABLE OF CONTENTS 


1.0 INTRODUCTION 

1.1 General 

1.1.1 Report Contents 

2.0 MAXIMUM LIKELIHOOD DEMODULATOR FOR THE SHUTTLE KU-BAND 
RETURN LINK 

2.1 Introduction 

2.2 Recommendations 

2.2.1 UQPSK Notation, Approach and Recommendations 

2.3 Derivation of Detection-Synchronization Schemes 

2.3.1 Broadband Channel 

2. 3. 1.1 Carrier Phase Recovery (Timing 
Assumed Known Perfectly) 

2. 3. 1.2 Timing Recovery (Carrier Reference 
Known Exactly) 

2. 3. 1.3 Non-Data-Aided Bit-Synchronizer 

2. 3. 1.4 Data-Aided Bit-Synchronizer 

2. 3. 1.5 Joint Recovery of Carrier Phase 
and Timing Information 

2.4 Synchronization of UQPSK Ku-Band Shuttle Return Link 

2.4. 1 General on Synchronization for UQPSK Signals 

- 2 Channel UQPSK 

2.4.2 3-Channel UQPSK Synchronization 

2.5 Conclusions 

2.6 References 


PAGE 

1 

1 

1 

6 

6 

; 

6 

6 

10 

10 

13 

24 

25 
33 

35 
' '38 

38 

42 

51 

53 




<=CinCc 


om 


TABLE OF CONTENTS (Cont'd) 


PAGE 


3.0 CARRIER SYNCHRONIZATION TECHNIQUES FOR UQPSK' SIGNALS USED 

ON SHUTTLE COMMUNICATION 55 

3.1 Introduction 55 

3.2 MAP Phase Estimator for Unbalanced QPSK Signals 56 

3.3 Unbalanced QPSK Carrier Reconstruction Loops Motivated 

by the MAP Theory for Known Bit Timing 60 

3.4 Performance Analysis of UQPSK Carrier Reconstruction 

Loops with Known Bit Timing 64 

3.4.1 Two Channel Costas Type Loop With Active Arm 

Filters without Hard-Limiters 64 

3.4.2 Performance of Two-Channel Costas Type Loop 

with Hard-Limiters 69 

3.4.3 Computation of the Loop-S-Curve and Its Slope 70 

3.4.4 Computation of the Noise Power 71 

3.5 Performance Comparison 77 

3.6 Unbalanced QPSK Carrier Reconstruction Loops Motivated 

by MAP Estimation Theory 82 

3.7 Performance Analysis of UQPSK Carrier Reconstruction 

Loops with Unknown Bit Timing 84 

3.7.1 Performance Analysis of a Two Channel Costas 

Loop with Passive Arm Filters 84 

3.7.2 Performance Analysis of a Two Channel Costas 
Loop with Passive Arm Filters and Hard- 

Limiters 90 

3. 7. 2.1 Loop S-Curve 91 

3. 7. 2. 2 Noise Spectral Density 93 

3.8 Performance Comparison 98 

3.9 References 102 

' 4.0 SHUTTLE LINK DEMODULATOR ANALYSIS OF SQUARING CIRCUIT 

MECHANIZATIONS USED IN SQUARING LOOPS 103 




■J!inC. 


om 


TABLE OF CONTENTS (Cont'd) 

4.1 Introduction 

4.2 Preliminaries 

4.3 Limiter Output Signal 

4.4 Squarer Output 

4.5 Weak Signal Suppression 

4.6 Effective Loop SNR and Squaring Losses 

4.7 Numerical Results 

4.8 Conclusions 

4.9 References 


PAGE 

103 

106 

no 

111 

112 

120 

125 

128 

129 


5.0 DETECTION OF SHUTTLE SIGNALS TRANSMITTED THROUGH THE 

TRACKING DATA RELAY SATELLITE (TORS) 130 

5.1 Introduction 130 

5.2 Communication Link Model 130 

5.3 Probability of Error for Shuttle' Signals Transmitted 

Through the TORS Satellite Repeater 135 

5.4 Confluent Hypergeometric Expansion Approach 140 

5.4.1 General Bandpass Nonlinear Channel 140 

5.4.2 Reduction to the Hard Limiter Case T43 

5.4.3 Asymptotic Values of P^ 144 

5.5 Backoff and Downlink SNR 147 

5.6 Error Function. Expansion Approach 150 

5.7 References 158 

6.0 ROCKET EXHAUST EFFECTS ON' SHUTTLE RF COMMUNICATIONS - 159 

159 
161 



6.1 Introduction 

6.2 Test Description 



TABLE OF CONTENTS (Cont’d) 


JUinC-t 


om 


PAGE 

6.3 Test Results 165 

6.3.1 Firing Results - Forward Antenna 165 

6.3.2 Simulation Results - Forward Antenna 165 

6.3.3 Model Results - Forward Antenna 165 

6.3.4 Firing Results - Aft Antenna . 165 

6.3.5 Simulation Results - Aft Antenna 165 

6.3.6 Model Results - Aft Antenna 165 

6.4 Analysis - Correlation of the Observation 184 

6.4.1 Firing Correlation - 184 

6.4.2 Firing Jitters 184 

Appendix 6A - Recorded Data 194 


7.0 


PERFORMANCE SENSITIVITY ANALYSIS TO DATA ASYMMETRY FOR 

UNCODED AND CODED SHUTTLE 8PSK KU-BAND RETURN LINK SIGNALS 197 

7.1 Introduction 197 

7.2 Data Asyimietry Degradation on the Linear White Gaussian 

Noise Channel 197 

7.2.1 Definition of Data Asymmetry 197 

7.2.2 SNR - Degradation 198 

7.2.3 Error Probability for Uncoded Signal 199 

7.2.4 Error Probability for Coded Signal 201 

7.3 Description of Bandlimited System 201 

7.4 Justification for the Detector 202 

7.5 Bit Error Bound 204 

7.6 Performance Degradation for Infinite Bandwidth 

Channel Case 209 

210 


210 



7.7 


References 



■<XinCc 


om 


1.0 INTRODUCTION 

1.1 General 


This is the final report on the SHUTTLE KU-BAND SIGNAL DESIGN 


STUDY performed for NASA Johnson Space Center under Contract No. 

NAS 9-14846 directed by William Teasdale. It represents the work 
accomplished during the period 15 May 1977 to 15 May 1978. 

The general objectives of the overall contract were threefold: 

(1) to perform a detailed analytical evaluation of the demodulation 
techniques for the Ku-band Shuttle return link; (2) to develop 
analysis and performance of uplink and downlink noises and nonlinear 
signal distortions on the detection of Shuttle signals transmitted 
through the TDRS satellite repeater; (3) to determine rocket exhaust 
effects on Shuttle communications. 

What follows is an overall description of areas addressed in the 
report, the results obtained and recommendations based upon the study. 

1.1.1 Report Contents 

This report is concerned with- addressing and documenting LinCom's 
findings on three main areas previously identified. The report is 
divided into 7 sections. 

Section 2.0 has addressed the problem of carrier synchronization 
and data demodulation of Unbalanced Quadri phase Shift Keyed (UQPSK) 
Shuttle communications signals by optimum and suboptimum methods. 

The Shuttle return link signal requires the implementation of a 
three-channel demodulator which is unconventional in the sense 
that the two channels (the lower data rate channels) are 
combined at the modulator by phase modulating the I and Q channel' 
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data streams onto a subcarrier. This subcarrier is then used to 
modulate the Q channel of a .carrier while the high rate I channel 
(50 Mbps) is phase modulated onto the I-channel of the carrier. 
Demodulation of this unconventional signal design is investigated 
herein using the theory of maximum likelihood estimation. Various 
physically realizable demodulator structures are derived under 
various assumptions regarding a priori knowledge of subcarrier and 
clock synchronization. Mechanization of these structures at the 
functional level are demonstrated and some generalizations are 
suggested'. The results are useful in defining the Ku-Band Shuttle 
demodulator for the, return link, i.e., the carrier, subcarner, 
clock and data detector functional diagrams. The ISI effect due to 
the bandlimiting of the channel is considered and recommendations 
to successfully combat this limitation are proposed. These 
recommendations can be useful when increased’ data throughput must 
be achieved while the quality of the performance of the system is 
to remain high. 

Section 3.0' addresses the problem of analyzing carrier 
reconstruction techniques for unbalanced QPSK signal formats. 

Certain signal designs for the Shuttle require carrier reconstruction 
from an UQPSK waveform in which the power ratio is variable between 
zero and six dB. It is also well known that a Costas or squaring 
loop cannot provide a coherent carrier reference for a balanced 
QPSK signal. Therefore, it is important to be able to predict 
the "optimum" carrier reconstruct! on technique for a varying power 
split between the inphase and quadrature channels of the signal and 
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compare the optimum method, with the Costas high rate tracking 
method. It is also of interest to compare this with that of 
conventional receivers for balanced QPSK, viz., the fourth power 
type, etc. Intuitively, some form of hybrid loop will arise as a 
consequence of the desire to track an unbalanced QPSK signal' as 
well as the more conventional QPSK signal in which the data rates and 
powers are balanced. 

There are various approaches to this problem^ including that 
of nonlinear filtering theory and the maximum a posteriori (MAP) 
approach. In what follows we shall derive the so-called carrier 
reconstruction technique based upon using the MAP approach. From 
the gradient of the MAP solution, closed loop' carrier reconstruction 
configurations are suggested. It will also be clear as to how the 
Costas loop fits into the problem of carrier reconstruction for an 
UQPSK signal. The theory is derived for two important cases; when 
the receiver has and does not have a priori knowledge of the clock. 

Section 4.0 is concerned wi-th evaluating the demodulation approach 
of the Ku-Band Shuttle return link for UQPSK when the I-Q channel 
power ratio is large,, i.e., 4:1. For this case the carrier can 
be recovered from the high power component only by the use of a 
squaring or Costas loop. Two different implementations of the 
carrier recovery loop for suppressed carrier BPSK signals are 
considered in this section. In particular, the tracking performance 
of the squaring loop in which the squaring circuit is mechanized- 
by a limiter-multiplier combination to produce the absolute value 
of the incoming signal, followed’ by appropriate bandpass-filtering. 
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is compared to that of usual squaring loop whose squaring circuit is 
the times-two multiplier which exhibits a square law characteristic. 

Section 5.0 considers the effects of uplink and downlink noises 
and nonlinear signal distortions on the detection of Shuttle 
signals transmitted through the TORS satellite repeater exhibiting 
both AM to AM and AM to PM characteristics. Probabilities of symbol 
errors as well as transition probabilities when the channel is 
considered as a discrete memoryless M-ary symmetric channel are 
derived as functions pf up/downlink signal-to-noise ratios and the 
characteristics of the TORS forward and return processors. The 
applications of the transition probabilities in the computation 
of bit error rate as well as the computational cutoff rate Rq are 
derived. The theory provided is needed in order to gain the 
necessary insight into Shuttle TORS simulations when the user 
constraints and other channel distortions are included into the model. 

Section 6.0 is concerned with determining the effects that 
Shuttle rocket motor plumes have on the- RF communications. LinCom 
has collected data taken at the Naval Weapons Center, China Lake, 
California-, at tests conducted by Lockheed Missiles and Space 
Company to determine rocket plume effects of the TRIDENT I on RF 
communication. Data was obtained on the effects of plume on 
amplitude and phase of the RF transmission at the GPS/SATRACK 
frequency of 1575 MHz. The three firing tests successfully 
measured attenuation, phase shift, attenuation jitter and phase 
jitter caused by the rocket motor exhaust. Empirical formulas 
are given for both the attenuation jitter and phase jitter in 



- 4 - 



terms of total RF attenuation. Knowledge of the effects of the 
TRIDENT I rocket plume on RF transmission will perhaps be useful in 
predicting the Solid Rocket Booster (SRB) plume effects on the 
Shuttle communication at S-band frequencies. 

Section 7.0 discusses the effect of data asymmetry on bit 
error probability which is of particular concern on the Shuttle 
Ku-band returnlink where large values of data asymmetry are 
anticipated. The performance sensitivity is analyzed for NRZ data, 
both coded and uncoded, and tabulated for various bit error 
probabilities. A theory which accounts for coding and channel 
filtering is presented and bit error bounds are derived for this 



-jCinCc 


om 


2.0 MAXIMUM LIKELIHOOD DEMODULATOR FOR THE SHUTTLE KU-BAND 
RETURN LINK 

2. 1 Introduction 

This section addresses the Ku-band service to Shuttle 
(return link) via mode one. In this mode three digital data 
sources serve to moduTate the return link Shuttle carrier 
(Fig. 2.1). Two of these three data sources are modulated 
onto a subcarrier to produce a UQPSK signal modulated subcarrier. 
This UQPSK subcarner then modulates the Q-channel of the carrier. 
The I-channel is modulated (high rate) directly with the other data 
source and these two signals (I and Q) are then combined to produce 
another UQPSK signal modulation for transmission to the ground via 
the TDRS. The effort on this task has been to identify optimum 
and suboptimum methods of carrier and subcarrier synchronization 
data demodulation techniuqes for this unconventional type of 
communication signal design. From this study a recommendation 
for implementation is made, 

2.2 Recommendations 

2.2.1 UQPSK Notation, Approach and Recommendations 

A simplifier UQPSK modulator is shown in Fig. 2.2. 
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The data source one is characterized by the modulation sequence 

{a^^} which can be real for BPSK modulation format {aj^=+l } or 

complex for a more general MPSK constelation {then aj^ is the 

complex root of unity or {aj^l = { = 0,1 . ,M-1}: 

The data rate is (T) ^ = R and the signaling waveform is p{t) which 

can be selected to comply best with the bandwidth constraints 

imposed by the channel. p(t) can be real or complex too (for 

SSB or VSB'type of modulation it is complex and real otherwise). 

/a =1 n = A 

The local oscillator waveform can be purely sinusoidal(^ -0 n 1 / 
or any other type of periodic waveform. The two signals s^(t) and 
S 2 (t) are combined to produce the signal s(t) which is 

s(t) = aj^Pl(t-kT.j)e'^'^0^| 

' k n 

+ /(l-a)P “nB|^P2(t-kT2)eJ'®6*| t^.^) 

' n k ' 

Here a is the percentage of the power P allocated to the channel. 

Using this notation to display all the relevant parameters of 
the modulation scheme, we arrive at Fig.2.1 which is the Shuttle 
orbiter transmitted functional model. The data on Ch. 1 can be 
MPT-L-jM-jS encoded or represent uncoded data, on Channel 2 it 
will be chosen to be either NRZ-L-jM.|S or Bi-<i>-L-jM^S and on 
Channel 3 it will be NRZ-L.jM.|S. The signal is observed in 
white Gaussian noise of single-sided power spectral density Nq 
watts/Hz. The satellite channel is assumed to be linear although 
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the existance of power amplifiers such as TWT in the link 
may cause the certain instances a concern about the channel 
linearity assumption . In addition, it is assumed that the 
channel is not bandlimited (which is a restriction) and this 
assumption can be easily removed by using a modification of the 
nonbandl imited solution as we shall later see. 

2-3 Derivation of Detection-Synchronization Schemes' 

2.3.1 Broadband Channel 

The signal at the output of the 3-channel modulator in Fig. 2,2 
(Mode 1 ) IS 


s(t) = ^ 


’^O.SP’ Re aj^p.| (t-kT.| )e 
* k 


I + Ii«|/ 07 r 5 F 


Re^ 


E 


n 1/ } 


>0^ 


+ /0.04P^ 


Im 


n t J 




( 2 - 2 ) 


In addition to a perturbing noise, the channel will add a random 
phase e(u) different for the carrier and the subcarrier and random 
delay t(u) different for each data stream. Thus the signal at the 
input of the demodulator will be 


y(t,0(u),T(u)) = s(t, g(u ),t(u)) + n(t,u) 


ORIQINAE PAGE IS 

QE eopjr; QUAUTja 


=v2 




+ /0. 1 6P’ Im 


- 1 n 1/ 


j(wgt+e.j (u) ) 
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+ /OTW Im 


^ n !/ 


j(no)s(.t+e2(u)) 


n k 

n(wQt+e^{u)) 


+ n(u,t) 


^-3) 


where the inclusion of (u) as an argument will indicate that the 
variable is random. 

Coherent demodulation of a carrier-modulated signal requires 
that a reference carrier with a precisely-controlled phase 
relative to that of the received signal be made available at the 
receiver. In the case of a double-sideband (DSB) signal, the 
effect of a relative phase error is only a loss in a signal -to- 
noise ratio in the demodulated signal. For other modulation 
formats, such as PSK and combined amplitude - phase modulation, 
which involve modulation of both the in-phase and quadrature 
components of the carrier, the problem is more severe in that a 
reference phase error will introduce co-channel interference into 
the demodulated signals (crosstalk). In such systems, it is often 
required that phase errors be held to less than a few degrees. 

Since the data signals are in the form of modulated synchronous 
stream (sequence) of pulses, it requires an accurate timing 
information for successful demodulation of the data sequence. 

The effect of an error in the sampling instants of the data 
detector outputs is an increase in the amount of intersymbol 
interference, which can drastically reduce the margin for errors 
when noise is present. The problem is especially severe when 
the data bandwidth is sharply limited to a value near the 
minimum (Nyquist) bandwidth for che data rate used. In this 
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situation, the timing error must be held to a small fracti 


on 


of the symbol interval (T) for satisfactory performance. 

In almost all situations, practical considerations require 
that carrier phase and timing information be extracted from the 
received signal itself. In some situations, it is practical to 
incorporate fixed periodic signals (such as pilot carriers or 
periodic data sequences, multiplexed with the information signal) 
in order to aid carrier phase and timing recovery. In other cases 
(as in ours) where the total power is sharply limited. It is 
important that all of the available power be allocated to the 
information bearing components of the signal. Thus, we restrict 
our attention here to the situation where no periodic signal 
components are available for phase and timing information. A 
unified approach to these problems is provided by considering the ML 
estimation of the unknown carrier phase and timing parameters. The 
ML estimator performance provides valuable bounds on the performance 
that may be attained with any other estimation schemes. Also, the 
structure of the ML estimator suggests a variety of suboptimum 
estimator structures, which may have practical implementation 
advantages. We examine two basically different approaches to 
phase and timing recovery problems by means of ML estimation. In 
one approach, the attempt is made to recover the unknown parameters 
without reference to any information about the data signal. In 
the other approach, generally called "data-aided" approach, we 
obtain superior performance by taking advantage of whatever is 
known at the receiver about the transmitted data signal. The 
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disadvantage in this approach (as we shall see) is that the phase 
and the timing estimators are coupled and that special means must 
be taken to start the scheme or decouple the phase and timing 
recovery. 


2. 3. 1.1 Carrier Phase Recovery (Timing Assumed Known Perfectly) 
Using the complex notation, introduced earlier, a general 

i 

carrier modulated signal immersed in noise can be expressed as 


or 


y(t) = s(t,8 ) + n(u,t) 

j(w«t-fe (u)) 

y(t) = Re{m(u,t)e } + n(u,t) £^4) 


where m(t,u) is the complex modulating random data signal and 
9(u) is the random phase introduced by the channel, and n(u,t) is 
assumed to be WGN(0,Ng/2) Ng/2 is the two-sided power noise spectral 
density. In complex notation n(u,t) = Re (c(t)+id(t))e^“0^ which 
is equi'-valent to the Rice representation of the noise in-terms of two 
low-pass processes c(t) and d(t). 

Restating the problem in complex notation will yield 


y(t) 


Re[{m(Lr,t)e-^^^'^^ +y (~u,\)} 


Y^(u,t) = c(u,t) + jd(u,t) 


(2,-5) 


The estimation of e(u) is based on observations over a time interval 
of Tg seconds. For effective performance, it must be assumed 
that e(u) remains relatively constant over this interval. The 
estimator performance will generally improve as Tg increases, but 
Tg must be held to a minimum in order to track more rapid phase 
fluctuations and to reduce the acquisition time. The likelihood 
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function in the case of WGN for the phase estimate e(u) will be 


given by 


L(e) = C exp 


' % 1 


[y(t)-s(t,e)rdt 


Since y(t) is independent of 9 , we can write 


L(e) = exp 


2N^ 


(m,m) 


’0 ^ 

In writing (2-7}the argument t and u were suppressed and the 
following defining relation for the inner product is used 


(z.y) ^ iRe{Y^,y^} 


= I-Re J 


Yz(t)Ty(t)dt 


■I 


' ( 2 - 6 ) 


(2-7) 


z(t)y(t) 


(2-8) 


(a) ML Estimator Not-Aided by the Data (Non~0ata Aided Estimator) 
In this approach the appropriate ML function to maximize 
(with respect to 0 ) is the expected value of (27) averaged over 
the random parameters in m(u,t). For a .BPSK-DSB with suppressed 
carrier, m(u,t) = +^1 at any time instant. Assuming that m(u,t) is 
a stationary Gaussian random process with autocorrelation K^(t-s) 
and since m(u,t) = ^ a^ij)^(t), the coefficients in the Karhunen- 
Loeve expansion of m(u,t) are independent Gaussian random variables 
and the expectation of L(e) is evaluated as follows 


m(u,t) = 2 a. 

i 


tcT, 


0 


(2-9a) 


> (2-9b) 
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Substituting (2-9a) into (2-7) yields 
Lq{0) = E{L(§}} 


Omginak Essa m] 
OF K)OK QMSSaa, 



exp 


2N 


0 





( 2 - 10 ) 


Evaluating the expected value, taking the logarithm and omitting 
the irrelevant constant terms yields 



E 


X ,«N„ 


(2-11) 


Noting that Aq (e) in (2-ll)‘can be written in the form of an inner 
product 

AqCs) = (6‘(jii>uj} (2'-12) 


the inner product of G-w and to where G-u means an operator G operating 
on oj(t} and 


G-(o(t) 


0)(t) = 


^ Xi+2No ^i 


Zl'^i'^i 7 
^ T 

from fe-10) 


(2-13) 


(2-U) 


The operator G is actually a system (in our case a linear system) 
related to the covariance operator K^(t) as follows: (see Fig.2.3) 
At moderately high signal -to-noise ratios, we can approximate 


G as a low pass filter with unity gain over the band where the 
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signal power spectral density is relatively large compared 
to 2Nq. Now if the signal process is sharply bandl united, 

(near the Nyquist bandwidth), we can say that G*to = tu (which 
means that the passage of nj(t) through a low-pass filter leaves it 
unchanged}^ then from (2-12) 

Ao{8) = / [Re(Yj,e''=®)]^dt (2-15) 


From (2-15) the structure of the receiver that evaluates Aq( 8) 
is evident. As shown in Fig.2.6, it involves a modulator followed 
by a squarer and an integrator. The closed-loop automatic tracking 
circuit which involves a voltage controlled oscillator (VCO) driven 
by an error signal proportional to 9i^/9e shown in Fig. 2.-5 is a 
direct result of the relation 



[Re(y e“^'®)][Re(-jy e‘-^'®)]dt 

y ^ 


Noting that 


(2-16) 


2Dle(Yye'j®)]CRe(-jYye'j“)] = Iiti{Yj,e'J®)^ (2-17) 

and that the integrator is actually a LPF, the scheme in Fig. 2. 5 
can be redrawn as in Fig.2.6. In Ffg.2.6, the two-lines mean 
complex signal and one line means that the signal is real. Since 

j(ajpt-f'0(u)) 

for BPSK, s(t) = Re{m(u,t)e } m = +1 at any t, an 

important interpretation of Fig. 2.6 can be made. 

To be able to track the phase we first must get rid of the 
modulation by squaring the signal and drive the VCO with the 
imaginary part of this operation. The obvious advantage of 
that interpretation is that it can be generalized for MPSK 
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signals by replacing the squarer with an M — power device and 
then we use the imaginary part of it as an error signal. As 
we shall see, this immediately yields the structure of the 
tracking device. For BPSK the structure of the actual tracker 
will be derived as follows. Denoting by a the real (inphase) 
channel and by b the quadrature channel 


(a+jb)^ = + j2ab - b^ (2-18) 

Im(a+jb)^ = 2ab (2-19) 

Thus, for BPSK, the driving signal for the VCO will be the 
product of the inphase and quadrature channel signals as shown 
in Fig. 2.7. This is the familiar Costas loop used to track suppressed 
carrier signals. For QPSK signals the complex tracker will be as 
shown in Fig. 2. 8. And the structure of the tracking device will then 
be derived based on the same method. The VCO will be driven by the 
imaginary part of the fourth power of the inphase and quadrature 
channel signals a and b respectively as follows 


(a+jb)^ = 

a^ t b^ + J4a\ - 6a^b^ - j4ab^ 

(2-20) 

Im(a+jb)^ 

= -4ab^ + 4a^b 

(2-21) 

Im(a+jb)^ 

= 4ab(a^-b^) 

(2-22) 


From (2-21 )-(2-22) two schemes for tracking QPSK signals can be 
proposed. The one based on (221 ) will be as in Fig2.9 and the 
loop based on (2-22) will thus be as in Fig. 2. 10. In Figs 2.9 and 2.10 
it is assumed that the gain factor of 4 can be incorporated in the 
gam of the loop filter. 
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Figure 2.8 QPSK Receiver (Complex Notation) 



s(t)+n(t) 


Figure 2.9 QPSK Receiver (Real Notation) 'OKIQINAIJ PAGE IS 

OFj ZOOK QUALrrXi 
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Figure 2.10 Alternative QPSK Receiver (Real Notation). 
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We consider now the case where the modulating signal m(u,t) 
is known at the receiver. For digital data signals, this can be 


accomplished by using the decision -directed approach whereby a 
reference data signal can be obtained during normal data trans- 
mission. This is known as the data-aided carrier recovery approach 
and uhe assumption is that the detected data symbols have sufficiently 
low error probabil ity so that the receivercan reconstruct an 
essentially perfect replica of the data signal m(u,t). 

The appropriate ML function is obtained directly from (2-7) by 
taking the logarithm and omitting the unnecessary constants, thus 
yielding 


Aq( 0) = ^ Re(Ty,me'^'®) 


0i\Q(e) 


1 


9§ 


2N. 


Re(v,,jme^®) 


(2-23) 


(2-24) 


■0 y 

Using the relations (228,2-23 and 2-24) can be readily written down 


J\3(§) = ^ Re[/yy(t)m*e'^k] (2-25) 

^ ^ Re[/yy(t)(-j)m*e“’^^dt] ( 2226 ) 

^ ' ORIGINAL PAGE IS 

For example, for a BPSK-DSB signal OF POOR QUALITY 


m(u,t) = ^ dj^p(t-kT) - +1; p(t) is real 

(2'27) 

Substituting (2-27) into (2-25) and (2-26) 'yields 
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(2-28) 




(9) = Re Y, - Re E 9|,-j-z(») (2-29) 


where 2 j^( 0 ) - J*Yy(t)p(t“kT)e ^"dt can be interpreted as samples 
Of a filter matched to p(t). 

Since SI'S to be obtained during the normal operation 
of the scheme, (2-29) can be implemented as shown in Fig. 2. 11. 

For the case of MPSK signals, {aj^} are complex constants 

\ = gj(2VM)mk m^ = 0,1,...,M-1 (2-30) 

Substituting ^-30) into ^-29), the MPSK decision-directed receiver 
will be 


'"o 

90 


(0) = Re Y. (-J)a^(e) = 


Re (-j)(a. +ja. )*(z. (§)+jz. (e)) 
k I Q *^0 

(2-31) 


( 

-^(0) = Re Y - 

S9 I u I Q 




(2-32) 


For a QPSK modulation { 3 ^^} = f+l,+]} = {a^ ,a^}. Implementing 
(2-29) then will lead to a receiver of the form shown in Fig. 2. 12. 

The delays in the circuits are included to compensate for the 
processing delays in the decision circuits (for a decision based on 
one symbol it will be one symbol time). The performance of ML receiver 


aided by the data offers no real advantage over the ML receiver not aided 
by the data for the case of BPSK, high S/N ratio and broadband channesl. 
However, for QPSK signaling and bandlimited channels, where additional 
degradations due to crosstalk and intersymbol interference become 



- 22 - 




- 23 - • 



















r£inCc 


om 


more severe than the degradations due to additive white noise, 
the data-aided ML receiver offers greatly improved performance 
This can he achieved by replacing the dotted box in Fig. 2.12 with 
a decision circuit, which bases its decisions on more than one 
sample-4"or example decision-feedback equalizer or Viterbi detector 
(VA), and replacing the delays in the loop to match the processing 
delays of the decision circuit. By properly designing the loop 
and using sophisticated decision circuits the impairments due 
to the bandlimiting of the channel (intersymbol interference) 
and those due to cochannel interference (crosstalk) can be 
effectively removed, thus the performance of the ML data-aided 
bandlimited QPSK receiver shown inFig2.12 asymptotically approaches 
the performance of the broadband non-data-aided QPSKML receivers shown 
in Figs. 2. 9 and 2.10. 

2'.'3.1.2 Timing Recovery (Carrier Reference Known Exactly ) 

The fact that we have the carrier reference ideally at the 

receiver allows us to convert down to baseband the modulated 

carrier waveform, leaving us with 

y(t) = x(t-T(u)) + n(u,t) (2-33) 

where t(u) is the random delay introduced by the channel, which 

IS to be estimated in determining the proper sampling times at the data- 

aided receivers. 


x(t) 

= a aj^P(t-kT) 
k 

(2-34) 

y(t) 

= m aj^p(t-kT-x(u)) + n(u,t) 
k 

(2-35) 


The likelihood function, derived in the same manner as for the 
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phase estimatTon problem is 


1{t) = exp J y(t)x(t-^)dt (2-36) 


2.3.1 .3 Non-Data-A1ded Bit-Synchronizer 

As before, in the derivation of (2-36) we. have assumed that 
Tq is sufficienly long, so that the random fluctuations in 
J X {t-x)dt term due to the data sequence can be neglected. 
a9 BPSK-DSB Signals 

For this case aj^ = +1 with equal probability, averaging over 
{ aj^} in L(^) yields 

Lq(^) = = E |c^ exp J y(t) a^p(t-kT-^)dt| 

Tq k 

= Ejl^ exp ^ J' y(t)aj^p(t-kT-^)dt| 


Aq(t) =ln(LQ(^)) 
a(t) = ^ In cosh 


1l,{t ) 


(2-37) 




(2-38) 


The evaluation of this quantity is performed by the structure 
shown in Fig,2-13. The close-loop tracking loop will seek 
the zero of the derivative 

= V tA -I (q,.(T))tanh q. (^)1 


(2-39) 


> 











And the actual implementation of (2-39) will be as shown in 


Fig 2‘.14. For the more general case when are assumed to be 
jointly Gaussian distributed with mean m and covariance matrix M 

= C y.-.y L(?,{a|^})exp |- i 5ZZ1 

_ CO ^ i j ^ 

Iq{x) = C expl - 1 ^ (2-40 

^ k j ’ 

The above integration is performed by completing the square on 
the exponent of the integrand to form a multivariate Gaussian p.d.f. 
(which integrates to one). The matrix Q is defined by 




+ 



r 


where 


r 


m 


I 


p(t-kr)p(t-mT)dt 


r((k-m)T) 


(2-41 


(2-4; 


u 


-2NQM'^ra 


(2-4; 


<1 

k 



y(t)p{t-kT-x )dt 


(2-4‘ 


The bias term u tend to zero for high S/N ratios (N^ 0) and is zero 

in any case for zero mean data. Since Q is a positive definite 
Toeplitz matrix, it can be factored (Q = v"^V) and thus we can 
write the log likelihood function as 


JjnCt 


om 


k 


Sk = ("-<» 

j 


(2-45) 


From(2-45)we see that the test statistic can be formed by filtering 

the sampled output sequence from the matched filter with a discrete 

filter {Vj^} to produce the output sequence {Sj^} or equivalently 

by modifying the matched filter so that its samples form the 

{Sj^} sequence- From the expression for Q, it is interesting to note 

that for high S/N ratio, the correlation in the data sequence (M 

matrix), has little influence in determining the ML estimator, while 

/ 

the intersymbol interference term (r matrix) has a major effect. 

The implementation of (2-45) is shown in Fig 2.15. Hence, the evaluation 

of Aq(t) is performed by the same structure as in Fig. 2.13 except 

that the In cosh(*) nonlinearity is replaced by (•) type of 

nonlinearity (which is the leading term in the series expansion of 
2 4 6 

In cosh ■*■••• )• ^0^ broadband channels and uncorrelated 

data, c((K-m)T) becomes zero for Kj^m and Q is diagonal, thus in 
forming the test statistic the digital filter is shifted out of 
Fig. 2.15. The closed loop bit synchronizer implements the equation 

k " k J Jo 

(2-46) 

And this can be pictorial ly represented in Fig. 2. 16. This scheme 
is the “narrowband" version of the "broadband" bit-synchronizer 
non-aided by the .data shown in Fig. 2,14. The difference is that a 
D.F. is included to combat the effects of the intersymbol interference 


)dt 
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and tanh (•) Is represented by (■) (which again is the leading term 
in its series expansion for small SNR). 


Instead of using M.F. for the bit-synchronizer, correlators 
can be used and this version of a bit-synchronizer will be as in 
Fig. 2.17. The phase of the VCC which controls the start and the 
termination of the integrate-and-dump circuits is changed every 
T sec by an amount proportional to the magnitude of and in 
a direction based on the sign of — ^ as computed during the 
previous KT seconds. 

As a further practical simplification, In cosh(-) can be 


approximated as follows 


In cosh X = 


xj >> 1 


-o- lx| « 1 


And thus Aq{t) can be approximated 




k °l\(. 


y{t )p(t-kT-T)dt| for large SNR 




y(t)p(t^lcT-r)dt I for small SNR 


T.(^) ^ (kT+^ < t < (k+l)T + 9 


Approximating the derivative — ^ by the difference function 


An(x + i At) - An(x - t At) 




















H >,0 


The practical approximation of the closed-loop bi t-synchromzer 
wi 1 1 be 


/ 1 yy 

NqAt ^ 


r 

f y(t)p(t-lcT-? + 4 )'** 

1 J'(t)p(t-k.T-T) - -^)dt - 


''t /'• Ax \ 
\U- -j) 


= < 


large SNR 


JL_ y 

V2NqAt k 


;/ 


y(t)p(t-kr-T- ^)dt 


to) 




\2 


y(t)p{t-kT-^+ ^)dt) 


Ail ~ / -f 

small SNR' 
(2-50) 


The closed-loop configuration using -^rSO) as error signal to drive the VCC 
will lead to the so-called early-late gate bit-synchronizer shown in 
Fig, 2, 17. For small SNR the absolute value nonlinearity is replaced by a 
squarer and the gam 77 ^ by — ^ and At is selected so as to optimize 

*^0 2Np 

the performance of the synchronizer. 


b. QPSK-D5B 

In a QPSK signaling system, one has available as input to the' 
synchronizer two independent BPSK signals. In particular, assuming 
perfect carrier synchronization, the output of the in-phase and 
quadrature carrier tracker of the QPSK receiver can be written, 
respectively as 


y,(t) = 

Xj(t,x) + hj(t) 

(2-51) 

yq(t) “ 

XQ(t,x) + nQ(t) 

(2-52) 

^ J (t , t) - 

y aj^p(t^<T-x) 

(2-53) 




It is a simple matter then to modify the BPSK bit-synchromzer to 
accomodate a QPSK input. Basically, each output of a carrier 
tracker (yj{t) and yg{t)) is processed by identical BPSK bit- 
synchronizers and the resultant error signals are added to drive 
the VCC as shown in Fig. 2.18. In practice the idealized bit- 
synchronizer can be replaced by more easily realizable early-late 
gate synchronizer with VCC again driven by the sum of the two 
I-Q error signals. 

2. 3.1. 4 Data-Aided Bit-Synchromzer 
a. BPSK-DSB 


For data-aided timing recovery we assume that {aj^} are known 
(from the decision at the receiver) and the relevant log-likelihood 
function is 



(2-54) 

(2-55) 


The implementation of (2-5^ is then shown in Fig. 2.19. 

The resemblance of this circuit to the closed-loop non-data-aided 
loop is striking. The only difference is in forming the error- 
correcting signal driving the V CC. In the latter case this signal 
is formed by multiplying the derivative of the convolved signal 
with the convolved signal itself and not with the estimated sequence 
(as in the data-aided case). This important difference will 
create ambiguity (false lock points similar to those in a carrier 
tracker or essentially the mean timing-wave frequency is higher 
than normal) whereas, the data-aided case, for with high 




-33c 




- 34 - 










probability and no ambiguities will exist. Thus, if we have a reliable 
estimate of the data sequence it can be used to greatly enhance 

the quality of the bit-synchronization. A circuit which actually 
generates locally the data sequence {q^} and then by monitoring the 
transitions in the successive data symbols produces a measure of 
the lack of synchronization between the input signal and the locally 
derived timing source. A block diagram of this loop (called data 
transition tracking loop DTTL) is given in Fig. 2.20 (BPSK-signals) 

This loop is fully explained in Ref. [2-1], pp. 442-467. A general ization 
for QPSK signal will again be two BPSK DTTL loops with an error signal 
driving the VCC which is the sum of the error signals of the two 
SPSK-DTTL loops, as shown in Fig 2.18. for Tow SNR, however, we can 
no longer claim that with high probability; thus, DTTL appears 

to be inferior to the absolute value early-late gate synchronizer 
shown inFig 2.17 (in case where {aj^} is uncoded sequence) 

joint Recovery of Carrier Phase and Timing Information 
It is a straightforward matter to extend the previous methods 
to determine simultaneously the ML estimates of t(u) and e(u) for a 
carrier modulated with a synchronous baseband digital signal. To 
illustrate some typical results we consider BPS'K-DSB case. In this 
case, the Re part of the signal Is 

y{t) = Re{m(u,t)e'^V“0^> + Re{y (u,t)e^"‘^0^} 

* 

".(t) - Ea^p(t-kT-x) 
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The log-likelihood function is 


a(tJ) = Re^Yy, ^ a^p(t-kT-T)e'^®j 

And the implementation of the closed-loop joint phase and timing 
tracking algorithm follows directly 


3? " ■ 2N„ S ^ p(t-kT-x)e'^'®) 

0 k 

If ^ ^ ^k Re(Yy,jp{t-kT-T)e^") 


The error signals (2-58)-{2-59) are used to control the VCC and the 
VCO, respectively, simultaneously. To evaluate the performance, we 
assume that with high probability (almost always) and make a 
series expansion of a( t, 8) around the points t=t, e=e. This leads 
to an approximate expression of the form 


T-T 

e-e 




where 


a = 


■Yr a,^P{t-kT-T)e^®. H a,, J p(t-kT-?)e^' 

0 \ k k 




2N7 X] ^k ^ 


b = 


-?w 


~ Ref^ a. p(t-kT-^)e'^*®,j V a. p(t-kT-f }e'^'®) 



inborn 


A 


B 


C 


0 \ k 9 t k 


1 


At. 


0 \ k k 


je 


2N. 


Rel 


f 

a,^p(t-kT-^)e'^® ,j XI 3k ^ p(t-kT-9)e 



(2-63) 

(2-64) 



(2-65) 


The size of the off-diagonal term represent the degree of 
coupling between the timing and phase estimates. However, since (2-65) 
actually represents the degree of interference due to the data, it 
can be eliminated by properly prefiltering the bandpass signal or 
equivalently by digitally filtering if the sampled recovered signal. 
Assuming that this is done, (2-65) can be made to vanish thus essentially 
uncoupling the joint estimation and performance of this algorithm then 
becomes as good as if either one of the parameters were known exactly 
previously. Things are now ready to apply the outlined theory for 
carrier and timing extraction to our specific problem-synchronization 
of the UQPSK Ku-band Shuttle return link. 

2.4 Synchronization of UQPSK Ku-Band Shuttle Return Link 
2.4'. 1 General on Synchronization for UQPSK Signals - 2 Channel UQPSK 
The simplified 2-channel UQPSK modulator shown in Fig. 2, 2 is 
characterized by the two in-phase and quadrature modulating 
sequences {aj^) and where each sequence has different data 

rate T^. , i = 1,2, and allocated power— that is, the I-channel 
modulating sequence { 3 ^} will carry a% of the total allocated power 
P and the Q-channel sequence {bj^} will carry (l-a)% of P. Therefore 
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the signal at the output of the UQPSK modulator at each time instant 
is 

j((«j.t+8(u}) 

y(t) = Re m(t,u)e ° {2-6i 

where 

m(t,u) =• yPJaj^p^{t-kT^} + j/PJ b^p 2 (t-kT 2 ) (2-6i 

|Pi(t)dt = E, ; fp 2 (t)dt = Ep [energy of the signals 

‘ sequences { aj^ and { b|^} ] 

The maximum-likelihood function to estimate the random phase e(u) 
will be (from (2-7)) 

L(e) = C^expRej^ (Y^(t)m(t)e'^®) - ^ (m(t),m(t)| 

(2-6 

where Yy(T) ~ ni(t)+Y^(t) is the complex envelope of the received 
signal, the complex envelope of the perturbing noise and 

the relations specified by (2-8)^ hold, i.e., 

^Re(Y 2 >Yy) = ^ Re y* Y 2 (t)Y*(t)dt = J z(t)y(t)dt = (z,y) 

T T 

Substituting (2-66a) into (2-67) yields 

L(;> = C, exp ReU- a.tj/PT b, )ej«) - i (m.rn)! 


(2-6 



L( 0 ) = C 2 exp Rej^ }exp Re ^ {y^j/FJb^eJ }j 


(2-69) 


Averaging over a. and b. in (2-69) (assuming that a. and b, are 


independent) yields 


^0 


• /PT 


Lq( 9) = cosh ^ Re{Yy^e'^®}cosh ^ Re(Y„,je'^®) (2-70) 


Taking the logarithm and omitting unnecessary constants will 
finally give 


A(e) = In cosh 




2N 


0 


[Re(y ,eJ® )] 


+ In cosh 


^ Re(Ty,oe'^®) 


(2^71) 


The closed-loop carrier tracking system will drive the VCO with 
the derivative of. a( 0), that is 


a\(e) _ 




2N. 


f RIQINAL’ page is 
POOR quality 


tanh ^ 510 ^ Re(Yy,e^®) 


’0 
+ tanh 


< . i5 

2[^ Re(y,^^ ) 


Thus 


0 




gg- Re(y^^^ ) 


90 


/P^ 

tanh ^ |Re(y^,,e^®)} 


/F^ 

2C 


’0 


/F" 


tanh ~ {Re(yyje^®)} 


/P~" 

2 


2Nq 


(2-72) 


From ^-72) the tracking loop implementation is derived and is 
shown in Fig. 2.21. One can interpret this scheme as follows; 
the loop processes the incoming I-Q signals by convolving it 
with filters ideally matched to Pn(t) and Po(t). Then it weighs 
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it by multiplying it with /P^ and /P^. Further nonlinear process- 
ing in the I-channel by tanh(*) and multiplication with the Q- 
channel signal wipes out the modulation in the I-channel and 
produces an error signal proportional to 6. Similar processing 
in the Q-channel by tanh (■) and multiplication with the I-channel 
signal is added to produce an error signal that will drive the 
VCO. For the limiting case when either P-j or ?2 zero, the 
loop reduces to the familiar Costas loop for BPSK signals. For 
QPSK signals (although there are many suboptimum schemes, 
for example, the loop in Fig. 2. 10) , it is clear that the UQPSK 
tracking loop can be modified easily for use with both QPSK 
and UQPSK by adding switches which will allow it to go from 
one configuration to another depending upon the modulation format. 

2.4.2 3-Channel UQPSK Synchronization 

The simplified UQPSK signal at the output of the modulator 

* 

(as shown in Fig 2. 1 mode'0 ) can be written as 

S(t) = Re ’(i/P^aj^Pj(t-kT^)+j[/P^b|^P 2 (t-kT 2 )-!^jVP^C(,i^(t-kT 3 )e ]e 

(2-73) 

The channel adds random phases 8, and' 02 to the carrier and 
subcarrier frequencies, respectively in addition to the WGN 
n(t), and random delays X 2 and x^. The received signal will 
therefore be of the form 


y(t) = Re{(/P^aj^p.|(t-kT.|-T^)+j’[/^I^P 2 (t-kT 2 -T 2 ) 


j(oj t+e, ) j(cDf,+ 0 -) 

+ J^C|^P3(t-kT,-x,)]e Me ° ^ )} 


3 ^3^ 


+ n(t) 


_(_2-74) 
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Assuming that timing is known at the receiver, the ML estimator for 
this signal will make use of 96 % of the total power P and' the theory 
for the 2-channel UQPSK developed earlier can be applied. 

Thus, estimating each random phase 9-j and 6^ by way of a 2-channel 
UQPSK demodulators, the proposed scheme for demodulating 3-channel 
UQPSK signal is shown in Fig. 2. 22. In the training period, the 
main loop will first lock on the carrier using the 80% of the power 
present in the {aj^} sequence. One way to ensure that the loop 
will lock on it is to send a training sequence known at the receiver. 
Next, the second loop will lock on the subcarrier and the heterodyned 
signal is converted down to baseband. Using then (aj^} and {bj^l 
as two channel UQPSK signal, the main carrier loop is aided by the 
additional power in {bj^} sequence to ensure proper tracking 
thereafter. 

There are 3 cases of interest, which will model the channel 
behavior accurately. Depending on each case, a different carrier/bit 
synchronization scheme which will perform best is proposed, 
a. Broadband Channel (No Intersymbol Interference) 

In this scheme, shown in the block diagram in Figure 2.27 
carrier tracking loops and the bit synchronizers are coupled and the 
bit timing was assumed known. This assumption can be relaxed although 
this will lead to a decrease in performance and therefore must 
be used only in the initial training period. The moment the carrier 
is recovered with sufficient accuracy, the bit synchronizer will 
be able to find the correct timing and aid the carrier tracking 
further, thus ensuring convergence of , the estimates of the phase 
and the timing to their correct values. The coupling between the 
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phase tracking and bit synchronization can be further relaxed 
by properly prefiltering the input signal (as shown in Section 2.3) 

In addition, there can be several variations of this scheme. 

For example, whenever we convolve the signal with a matched filter 
we can use integrate and dump filters and vice versa. Also, 
instead of early-late gate bit synchronizer, for example, the 
DTTL synchronizer shown in Fig.2.20 can be used. Also, if 
more reliable data detection is desired, the data detection 
boxes can be -replaced by more sophisticated decision devices 
decision feedback equalizers or ML v^equence estimation imple- 
mentation by the Viterbi algorithm . Little gain will be 
achieved, however, by doing this replacement for our case — 
broadband channel, high SMR, where the performance of the scheme 
in Fig2.27 is nearly optimum. For the case, however, when there 
is intersymbol interference due to bandlimiting in the channel, 
a significant improvement in the performance can be achieved by 
basing the implementation of the block diagram in Fig. 2.27 on the 
data-aided approach as shown in Fig, 2. 28. 
b. Moderately Bandlimited Channel (Small IS I) 

I As mentioned before, to ensure that there will, be no 

drop in performance due to ISIthe carrier tracking loops should 
' / 

be data-aided. This will ensure best possible carrier 
reconstruction and will keep the crosstalk and distortions as 
low as possible. The bit synchronizers can be left non-data 
aided -Si nee they will operate well enough in a small IS I environment, 
yet keeping the coupling with the carrier-reconstruction loops to 
a minimum. This, in fact, will allow a, fast start-up and 
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convergence of the whole synchronization system since no special 
circuits for a start-up will be needed. Although in general for 
low SNR the data-aided bit-sync is superior, in the particular 
case of small ISI their performance will be very nearly equivalent 
and thus the increased complexity of providing start-up circuits 
does not justify the use of a data-aided bit-synchronizer for 
this case. 

To find the structure of the data-aided carrier tracking 
loop for 3-channel UQPSK we first derive it for the 3-channel 
UQPSK as in Sec. 2. 411 and then generalize it for the 3-channel 
-■'^ase. The ML function for the 2-channel case is derived in eq.(2-69) 
^ However, since we are trying to come up with a data-aided scheme, we 

‘ need not average over {aj^> and {bj^} (in fact, they will be 

assumed known or estimated). Taking the log in" (2-69) yields 

A A 

. " I’® i ^ (2-75) 

The VCO will be driven by the derivative of A(e) with respect to 
e, namely 


aie) - Rp 
30 2N, 


+Re 5 ^ (yy.-^(-bfe)e-'®) 


1 




0 


(2-76) 

The implementation of (2-7^ will lead to the circuit in Fig. 2.23 
which is the two-channel UQPSK data-aided carrier reconstruction 
loop. 

Following the considerations which led us to the 3-channel 
ML non-data aided UQPSK in Sec. 2.4.1 the 3-channel ML data-aided 
carrier/tracking loop is shown in Fig.2.24' and the complete 
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carrier/bit synchronizer is pictured in Fig.2>-28', 

t 

c. Severely Bandlimited Channels (Large ISI) 

When the channel, however, is severely bandlimited, to ensure 
distortion-free and reliable communication, we must employ the 
data-aided approach to both carrier tracking and bit synchronizing 
circuits. This will lead to the circuit in Fig.2.29. In this 
block diagram, the bit-synch circuit will be data aided, thus 
the DTTL circuit fromFig 2.20 is appropriate. To ensure best 
possible operation, the simple decision circuit (quantizer) in 
this scheme can be replaced by a VA box or use the already 
estimated sequence {aj^} from the carrier circuit to form the 
transition signal Since a delay is required to produce 
an estimate {a^^} in the VA box, the corresponding delay must be 
inserted in the lower branch of the bit-synchronizer leading to 
a circuit shown in Fig 2. 25. Now the derivation of the bit-synch 
depends on which are derived from the carrier tracking loop, 

f. 

thus tne estimates (e,t) are coupled. Uncoupling can be 
- . accomplished by sending a training sequence. or filtering the 
input of the bit synch properly. If this is to be avoided, 
however, then the circuit in Fig. 2.25 can be modified by using 
a separate VA box to replace the simple decision circuit. Since 
the decision is now based on several data bits, an appropriate 
delay must be inserted in the lower arm to compensate for the 
processing delay. Thus, using either one of the schemes in Fig. 2.25 
or Fig»2.26can serve as an appropriate data-aided bit-synchronizer* 
box for the strongly bandlimited channel in Fig.2.29, 
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Figure 2.25 Open Loop Bit Synchronizer for a Severely Bandlimited 
Channel . 
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2,5 Conclusions 

A theory for carrier and bit-synchronization for QPSK and UQPSK 
is presented. From the theory, circuits performing carrier track- 
ing and bit synchronization are derived. Three cases of interest 
are examined: broadband, moderately narrowband and severely 

bandlimited channel cases. Recommendations considering the 
tradeoff between the simplicity of the implementation and the 
reliability of performance are given. Specifically, 3 synchronization 
schemes are proposed: for the broadband channel a synchronization 

scheme using non-data aided carrier/bit synchronizer as shown 
in Fig, 2 -27 is recommended. For the moderately bandlimited 
channel a carrier reconstruction loop will be data-aided while 
as bit-synchronizer a simpler non-data aided circuit can be used, 
as shown in Fig.2,28. For the severely bandlimited channel, however, 
both the carrier tracking loop and the bit-synchronizers must be 
data-aided, thus leading to the block diagram in Fig.2.29. This 
scheme will eliminate the impairment imposed by the bandlimiting 
of the channel and will ensure high performance even in most 
difficult conditions, once started. 
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3.0 CARRIER SYNCHRONIZATION TECHNIQUES FOR UQPSK SIGNALS 

USED ON SHUTTLE COMMUNICATION 

3.1 Introduction 

An unbalanced QPSK (UQPSK) signal format is a QPSK signal 
where the two quadrature components of the carrier have unequal 
power and are modulated with independent binary data streams of 
different rates. This transmission scheme is most useful when 
the data of two different sources with possibly varying rates have 
to be transmitted in the same frequency band. This signal design 
will be used in the development of the Tracking and Data Relay 
Satellite System (TDRSS), Shuttle Communications, and the Navigation 
Satellite Timing and Ranging Global Positioning System (NAVSTAR GPS). 

Carrier synchronization techniques for UQPSK signals in which 
the power split in the in-phase and quadrature (I-Q) channels is 
variable and the I-Q arm data rates are also variable and therefore 
of interest. In fact, certain signal design concepts for the TDRSS 
require carrier reconstruction from an unbalanced QPSK waveform in 
which the power ratio imbalance, say X, is variable between zero and 
six dB. It is also well known that a Costas or squaring loop 
cannot provide a coherent carrier reference for a balanced QPSK 
signal. Therefore, it is important to be able to predict the 
"optimum" carrier reconstruction technique for a varying power 
split and data rates between the inphase and quadrature channels 
of the signal and compare the optimum method with the performance 
of a Costas loop tracking the high power, high data rate channel 
only. It is also of interest to compare this with that of 
conventional receivers for balanced QPSK, viz., the fourth power type, 
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etc. Intuitively, some form of hybrid loop will arise as a consequence 
of the desire to track an unbalanced QPSK signal as well as the more 
conventional QPSK signal in which the data rates and powers are 
balanced. 

There are various approaches to this problem including that of 
nonlinear filtering theory and the maximum a posteriori (MAP) approach. 
The carrier reconstruction technique is derived based on the MAP 
approach outlined in Ref. 3-1. From the gradient of the MAP solution, 
closed loop carrier reconstruction configurations are suggested. 

It will also be clear as to how the Costas loop fits into the 
problem of carrier reconstruction for a UQPSK signal. 

The theory is derived fro two important cases: In case one the 

reel ever has no knowledge of the clock whereas in case two the 
receiver is assumed to have knowledge of the clock. This latter 
case is not too restrictive an assumption since high quality 
frequency control can be used (as e.g. for the TDRSS) so that the 
carrier frequency uncertainty is small when compared with the link 
data rates. Thus the beat note, which serves to lock the carrier 
loop, is of sufficient quality to lock the bit synchronizers 
that provide the clock prior to carrier lock. 

3.2 MAP Phase Estimator for Unbalanced QPSK Signals 

An unbalanced QPSK (UQPSK) signal is defined as 

s(t) = d.j(t) cos(wQt+e) + ^2P^ sin(topt+e) ^ (3-1) 
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where P-j and ?2 represent the power in the two quadrature components, 
d^(t) and d 2 (t) are two independent binary signal waveforms, Ug/Zstis the 
carrier frequency and 9 is the carrier phase. The problem is to find the 
maximum a posteriori probability (MAP) estimate of the phase 6 when s(t) is 
observed over a time interval T and is immersed in white Gaussian noise 
of two-sided spectral density Ng/2 


x(t) = s(t) + n(t) 0<t<T 

The derivation follows the same outline as for a BPSK signal (Ref 3-1 
p. 524-530} and will therefore be kept short. 

We assume that data source 1 (generating d-j(t)) is clocked at times 
kT^ ,k=0-,l ,2, , . . ,Kj and data source 2 is clocked at times kT 2 > k=0,l,2,..., 

K 2 * Let ,k=l ,2, . . . ,K^+K 2 be the ordered set of times t and 5 . Let 
tg = 0 be the beginning and tj^ +k +] " ^ observation in- 
terval [0,Tj. The data in (3-1) are then constant in any interval »^j^) • 

The observed signal x can be expanded into orthogonal base functions cos upt, 
sin ugt in each such interval: 

x(t) = aj^cos ugt + bj^sin u^t + nQ(t); tj^_-j<t£t|^ (3-3) 


where ng(t) is an additive noise process which is uncorrelated with the signal 
and the noise components contained in and and can therefore be 
neglected for the estimation. Conditioned on 9 ,d.|(t) and d 2 (t) the 
random variables are independent and Gaussian with mean 

£{ aj^i .dg.e )]’4iP-]d^ (t) cos e+'iiP 2 d 2 (t)sin 8; t (3-4) 








E{ e+'l2P2d2(t)cose ;t e{t^_^ ,t,^3 (3-5) 

and variances where<itj,=t|^-tj^_-, . 

For brevity in notation we define the vector 


and let p(e|c ) be the a posteriori probability of 9 so that Bayes 
theorem allows us to write 


(3-6) 


P(e) 

P(6|c) = p(^ 


[d^,d2} 


p(c i d^(t),d 2 (t)^ ) dd^dd 2 


(3-7) 


For a uniform probability density of 9 the HAP estimate of 9 is the value 
which maximizes the integral above. The integrand is easily written as 

K.J+K2+1 At|^ j . 

p(c!d^(t).d2(t),9) = C({tj^»exp IT ^0 ' 


exp 


K^+<2+1 Atj^ 


k=l 


k=l 

K^+K2+1 it 

Z 

k=l 


a^»j2P^di(tk)cos e- N„ “j 


K,+K,+l « k,+K2+i ^ 

E E N. ‘>k^2-‘2'VPPPP 


K^+K2+1 it,^ 


(3-8) 


where C({tj^}) is a constant independent of 9. Collecting terms in the sums which 
contain the same data symbol d-j, d 2 and then averaging over the data yields* 


kT, 


p(£|e) = C 2 cosh! Nq \f2pT| ( x(t)cos(aiQt + e)dt 

(k-l)T, 


K 


r cosh I Nq \[2P^ x(t)sin{ WQt+9)dt 

•- (k-l)T2 


kT, 


^C,f(9) (3-9) 


*To simplify the notation we assume T = K-jT^ = *^2^2' 
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where cosh x is the hyperbolic cosine function. The optimum estimate for 
0 maximizes f(e) or equivalently In f(e), i.e., 


In f{e) 


^1 

E 

k»l 


kT, 


In cosh 




S. 


x(t)cos( 03 gt+ 0 )dt 


(k-l)T, 


+ 


X) In cosh 
k=l 


J- 'Jzpl \ x(t)s1n( a)nt-»6 )dt 
LO ^ (k-l)T2 ° 


(3-10) 


We can interpret (3-10) as follows. For each value of 0 , one correlates the 
received signal x(t) with the stored local replicas ( 2 /Ng)'T^ •cos(ugt+ 8 ) 
and (2 /Nq) \ 2P-^ sin(<^gt+ 9) over respective intervals [(k-1 )T^ ,kT^],[(k-l )T 2 ,kT 2 ] 
The result of the correlation in the respective sums of ■ (3-10)- are 

added and weighted in accordance with the log hyperbolic cosine. These 
random variables are accumulated over K-j or K 2 successive symbol intervals 
and stored. The phase 9 yielding the largest accumulated value is 

A 

then declared as the MAP estimate, say 9. Unfortunately, 0 takes on a 
continuum of values between -^and and hence to apply the above tech- 
nique in practice for eachTsecond observation interval one must quantize 
©into a number of levels compatible with the size of buffer storage that 
one can economically provide and the accuracy of the estimate required. 
Furthermore, 9 is usually a slowly varying process when compared with the 
center frequency uq so that some form of updating or tracking the changes 
in 9 is desired. To accomplish this we now propose closed loop systems 
which can be motivated by the MAP theory. 
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3.3 Unbalanced QPSK Carrier Reconstruction Loops Motivated by the MAP 
Theory for Known Bit Timing 

Any successful carrier reconstruction technique 
must provide for acquisition as well as tracking of the slowly varying 
(now assumed) process 8{t). Thus a closed loop concept which employs 
the phase-locked loop principle is of interest. It is possible to pro- 
pose a closed-loop system, incorporating the development given in Sections. 2 
that would continuously update the estimate of the input phase e • This 
can be accomplished by differentiating f( e), or equivalently In f(e), 
with respect to eand equate the result to zero. The solution of this 
equation represents the peak of the a posteriori probability density func- 
tion and hence the MAP estimate. Hence, 


g( 0 ) =t^Onf(e)) 

- E 

k=l 


Kj 

E 

k=l 


tanh 


6=6 

kT 




1 


L^O 


x(t)cos((Unt+9)dt 


tanh 






{k-l)T 
kT. 


J *^0 


1 ^) 


,kT, 


1 


(k-l)T 


x(t)sin(uQt+e)dt 


1 


2 j x(t)sin(o!Q+e)dt 
(k-l)T2 



.kT. 


J Nq ^ 


A 

x(t) cos(uQt+e)dt 


(k-l)T, 


(3-11) 


By using g{e) as an error signal in a tracking loop to control the phase 
of a variable phase oscillator one arrives at the closed loop configuration 
in Fig.3.1_. Basically, the receiver takes the in-phase and quadrature pro- 
jections of the incoming signal components and performs correlation inte- 
grations over the respective bit intervals. The voltages are held for 
time intervals equal to the symbol durations of the respective data 
rates and weighted in accordance with the power split in the channels. 
Further processing includes the hyperbolic tangent weighting and then 
accumulating the samples over and K 2 symbols in each channel. The 
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difference of the two arms is taken and this provides for a phase esti- 
mate which is used to control the phase of the local oscillator. 

Two points are worth noting here. First, the theory has suggested the 
optimal weighting, for a given power split at the transmitter, of the 
decision voltages in accordance with the square root of the signal power 
to noise spectral density and not "P/Nq" itselfl This fact is important 
when the effective loop gain, which sets the loop bandwidth, damping and 
static phase error becomes of concern in the design. Furthermore, when 
P-j equals zero, I.e. , a binary phase-shift keyed signal is transmitted, 
the upper arm of the loop opens and the lower loop reduces to a configura- 
tion reminiscent of a Costas loop. The same comment is true about the 
lower arm when P2=0. On the other hand, when P-j=P 2 ^nd T-j = T£ we have 
a configuration which will reconstruct the carrier for purposes of de- 
modulating QPSK. 

Using analog circuitry the loop is probably difficult to implement 
owing to the requirement of implementing the nonlinearity tanh x; how- 
ever, it is interesting to examine the carrier recontruction loop which 

results when the signal -to-noise ratio is assumed to be high and low. 

The hyperbolic tangent of x can then be approximated by sgn(x) or x, respectively. 
Using these approximations in (3- IT) leads to the implementation of Fig. 3.2 

where G{x) = sgn(x) for high SNR and G{x) = x for low SNR. Note that the 
number of demodulators has been reduced to two and that the digital filters 
(accumulators) have been replaced by an analog filter in this implementation. 

These loops will be referred to as two channel Costas type loop with active 
arm filters. Notice that since the accumulators required by the MAP approach 
have been removed here the optimum values of the channel weighting constants are 
K-, = 4 P^/(Nq T^); K 2 = 4P2/{N^T2) (3-12) 

where T^.= RT^ is the bit length of d^(t). Equivalently, the optimum ratio 
of the two gains is 


K 1 /K 2 = (P^T2)/(P2T^) 
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3.4 Performance Analysis of UQPSK Carrier Recontruction Loops with 
Known Bit Timing 

3.4.1 Two Channel Costas Type Loop With Active Arm Filters Without Hard-Limiters 
Consider the loopof Fig3.2 with G(x)=x and the input signal x(t) = s(t}+n(t) 
where s(t) was defined in (3-1) and 


n{t) =JT [N^(t)cos{(OQt+9 )-Ng(t)sin((i3Qt+e)] 


(3-14) 


Since 6 is constant, or at least slowly varying, the correlation between 
and can be neglected. They are assumed to be baseband Gaussian 
random processes with a flat spectral density NqW/Hz (single-sided) 
in their bandwidth . Neglecting double frequency terms the output of 
the phase detectors (gain'll ) becomes 

£^it) x(t) \12 Kq COs(a)Qt+0) 


= C^d^ (t)+N^(t)]cosvKt) + [jP^d 2 (t)-N 5 (t)]sinv(t) 

x(t)4HQ sin(o)Qt+e) 


(3-15) 


IKqK^ [P^tl^(t)+N^{t)] Sinv>(t) - [4p‘2d2(t)-Ng(t)]cosv>(t) 


(3-16) 


where v^e-e is the loop phase error. If the phase error process is 
slow we can assume that<pstays constant over one integration time 
(for both transmission rates). Setting KqK^^ = 1 this yields, see Fig.3..2 

Zci(t) = T. 'Jp"^d^.(t) cosv»+ T- \[P^ d2- (t) sinv’ 


+ N^^.(t) cosv?-N^^(t) sinv? 

Zsi(t) = T. d-]^*(t ) sin^-T^ d 2 .j (t)cos 

+Nc.j(t ) sin v»+N^. (t) cos 


(3-17) 


(3-18) 
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where 




kT. <t<(k+l)T^ (3-19) 

kT^ <t< (k+l)T^ '(3-20a) 

kT^<t< (k+l)T^. (3-20b)‘ 

' 2 
^ci’ ^si’ constant functions with variance a = NqT^/2. 

Since N ,N are band! imited, successive outputs N ,N of the integrate-and-dump 
c s ^i ^i 

circuit are correlated. However, if the noise bandwidth is considerably larger 
than the faster signal transmission rate this correlation may be neglected. 

This will be done here. Note, however, that may be correlated with^ 
several successive samples of vice versa due to the different 

integration times for the two channels. Since the quadrature components 
are independent we can write for the crosscorrelations 


£1 


^ J' 

' / 1 •'n \ -» 


d^(s)ds. 


(k-l)T. 


kT, 




N^{s)ds, 


(k-l)T^ 

kT, 


"^(Ods, 


rif-1 IT. 


E{N^j^Ngj}= 0, j,k = l,2 

E{Nci(t)Nj,2{s) 1= E{N5^(t)Ng2^s) 


(3-21) 

(3-22) 


The outputs of the multipliers are easily shown to be (see Fig^3.;2 
2,(t) = K,Z^,(t)Z3^{t) 

=■ j{T2[p^d^.(t)-P2d^,{t)]+2T.[,fF^d,,UXj.(t) + F2d2i(t) Nj,(t) ] 

+ N^.(t ) -N^,.(t) )■?. si'n(2»>) 

T?JP7Jd„(t) dj^.(t) + T.[iF,d,,(t) Nj^(t) - Pj dj.Ct) N^,(t) ] 

+ N^.(t)Nj,.j(t)} K. cos(2^) (3-23) 
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The loop phase estimate of the incoming phase is derived from z-j and Z 2 '- 
- K/(P) 

~P~[z,(t)-22(t)] (3-24) 

♦ 

For 5=0, the case of greatest interest here, and Kj^=KyKj^ we find 
2pv> + (K^T^P^ + K2T^P2)F(p)sin(2v’ ) 

- F(p)jE[K^T^P 2 d^t) ]sin{2v>)}= F(p)N(t) (3-25) 

where N(t)= N.|(t)“N 2 (t) and 

N^(t) = N.^(t) + N.2(t) + N^3(t) i=l,2 (3-26) 

N^-T(t) = Ki[N^(t) -N^.(t )]sin(29)+ 2K.N3^.{t )N^^(t )codZ<p) 

\ / a ^ 

N- 2 (t) = 2K.T^[>Jp^ d-j^. (t) Nj,.(t) +^IP^ d 2 ^ (t) N^^.(t ) ]sin(2v>) 

+ 2K.T.[Jp^ d^.(t) Ng.(t ) -f^d 2 ^(t), N^^.(t) ]cos(2v?) (3-27 Jd) 


Nj3(t) - JvT<l„(t)d2,(t)cos(2p) 

{3-27c)‘ 


We note that the three terms m (3-27) are zero mean and uncorrelated for 
a given V5. This simplifies the computation of the spectrum needed in 
obtaining loop performance in- noise. 

Defining the signal component 


S'(0)=^ 


K,T?P,+KoToP2-K^T^P2E{d2^(t^K2T2P-,E{d^2(^)^ 


'Vr l‘'2‘2^ 


(3-28) 


and linearizing the equation of operation the phase error is given by 


^ H ^(p) [N{t)/S'{0)l 

where 

C S'{0)]F(s) 

H ^(s) = 


(3 -.29) 


(3-30) 


s -H [S'(0)]F(s) 





Assuming that the bandwidth of N(t) is much wider than the bandwidth of 


H <^{s) the variance of ^ can be approximated by 
) 


.(3--31 ) 


where is the one-sided loop bandwidth and is the one-sided spectral 
density of the equivalent noise N(t). 

As noted earlier, "ck (t), Ngj^{t), k=l,2 are piecewise constant with 
discontinuities at the signal clock times. Hence if the clock phases are known a 

N 2 are nonstationary. For the nonstationary correlation function we 
find [3-3] 


» 4Kf<7^Ca2+T^P^+T^P2E(d^^]] + T^PTp2ECd^^]} A^(t^,t2) 
+ ^ A2(t^,t2) 


2r22 


-4 J <1 ^ ^2 > ^1 ^ ^1 ‘'”2 ^ d , T ^ ^ 4 ’ ^1 d,., ^^2 ’ ^1 ^ ^ c 2 ( ^ 2 » ^1 ^ 


’cr'c2 " ' ' ^ ' ^iri2 ^ ‘ *• '^2r22 

r2^2. 


+Ti 12?^ P2Rd.^ d^ 2 ’ ^1 ^ ^d2-| d2j ^^2 * ^1 ^ ^ 


(3-32) 


where 


ORKawAi,- pa™ „ 

OF pnon IS 

^ qvauty ; 


A.j (t.j *t2) ^ 


1 if t-j and t 2 in the same interval (kT^ ,(k+l )T^ ) 


0 otherwise 


■^S-33) 


To keep the computation simple we assume that T 2 =nT-j where n is an integer 
and that the clock of data source 2 is derived from the clock of data source 
1, i.e. each clock pulse of d 2 coincides with a clockpulse of d.j. (This can 
be shown to be a worst case assumption. For a more general result see [3-3].} 

The variances are then found to be • Computing the 

correlation functions is straightforward but tedious [3-3], The result is 
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1 Oct+TOiT-j 


R. . (t,t+T) =\1 t>T, , nT,<t+T<2nT, 

^2r22 f 

to otherwise 




(t.t+x) 


(■3--34) 


Sl“l2 “ " ''''2l‘‘22 


(t,t+T) 


Substituting these functions into (3-32) for t.j=t, t 2 =t+r and averaging 
over t yields the stationary correlation function which can then be 
integrated over t to yield the spectra! density of the noise at ( 0=0 


*S^{0) = N„T^{K^[Ng + 2T,(P,tP2)3 +l<2n^C''V^''''', (P^+nPj) 

- 2nK^K2[NQ+2T^(P^+nP2)]| 

-f4T^P^p2[)(^ + - ZK^KgO^j (3-3,5) 


The phase error variance is found from (3-31) 


No 


2^3 


■V ■ "o\ 




/ 2PtPoT,\ o . 

ZK^Kgn]^- (P2->- --47/ 


[K^P^ +K2P2n"^ - K^P2' nK2P^3' 


(3-36) 


We note that the performance of a single channel Costas loop can be 
found from this result by setting either K-j or to zero. To find the 
tracking performance of a Costas loop for a BPSK signal P 2 and K£ can be 
set to zero. The result then becomes independent of n. 

The squaring loss is defined as 


- 68 - 




2 

<r 


Vi 


^(P,*P2) 


- [S'(0)] /P 
"o''"o 


X Co 


m 


(3-37) 


A ^ ‘ 

Introducing P = H 

constants can be rewritten as 


and using the optimum ratio of the 



[l-2X+nX^]^{1+X)“^ 



(l-2x+nX^)+X 





(3-38:) 


This loop IS unable to track Da lanced QPSK since for n=l 






(3-39) 


which is zero for X=1 

3.4.2 Performance of Two-Channel Costas Type Loop with Hard-Limiters 
This implementation is represented by .Fig. 3. 2 when G(x) = 

sgn(;^. Due to the nonl inearity in the loop a slightly different 
approach must be used here to find the performance. 

The dynamic phase error z{t) can be split into a random and a non- 

random part for a given phase error ^ 

z(t) = E{z(t)|P} + N(t,<P)=S(«F)+N(t, «P) (3-40) 

where S{<p) is the loop S-curve and N(t,vJ) is a zero mean white noise process. The 
loop equation can then be written as 

_JL |S(¥*) + N(t,^)} (3-41) 
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The loop equation can be linearized by setting 

S(v?) w S(0} + S' (0)^ 

N(t,v?) « N(t,0) 


JZlnC 
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(3-42) 

{3-43:)'- 


where S(0) =0 and S'{0)= 


jv>=o • rewritten as 


S'(0)K/(p) N(t,0) ^ ^ N(t,0) 

p +S'{0)KyFM “S'W S'(0) 


(3-44).. 


which is identical in form with (3-29) , hence the phase error variance 


follows from (3-31). 

4 . ■^ Computation of the Loop S-curve and Its Slope 
The loop S-curve is defined by 


SM = \|i^S^(<^) -^S^[ ¥>) ^ (3-45) 

where 

S^{<p) = E{z^^(t)sgn[z^^(t)]|‘p} (3-46)' 

= Ejz^2(^)s9"tZs2(^^^l'^} (3-47) 


From (3-17,3-18) it is clear that z^pZ^I are independent Gaussian random 
variables when conditioned on d-j-jjdg-j. Hence the conditional expectation 
of (347) can easily be found 


S^(v’) = E (Zg^sgn(z^^)} = 


T !('rr- sinv?- \|p 1 cos¥’)erfi\fK cos¥»+ sinv> )— 




+ ( s1nv» + ^ cos<p)erf 


(>TF| cosv’ sinv») — 




(3-48) 


1 J 


and similarly for the other channel 
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= E{2^2Sgn(Zj2>’ * ^ 


n -j 

I] jh(r -7 cosv= + ^ sln^J 

k=0 '• 




COS^P - \fp^ 


(3-49) 


S =>[^ Si{<p) ^2^^^ is plotted fn Fig-3.. 3 for n “1,2. It must be 
pointed out that this loop implementation has multiple false lock points 
and that the slope of the S-curve depends on the signal -to-nolse ratio. 

Using (3-48), (3-49) ,'eq. (3-45) can be differentiated to find the slope 
S' at zero 



(3-50) 

where we have used the optimum ratio of the gains (3-13). 

3.4.4 Computation of the Noise Power 

Under the assumption of a small phase orror the signals in the four 

arms can' be approximated by 



(3-51) 

Zs,(t) • -T, .[P7<i2,(t) +Nj,(t) 

(3-5Z) 


om 
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Figure 3-3a, Typical Loop S-Curves for Balanced QPSK. 
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In these terms the dynamic error signal can be written as 


X/nC 


om 


z(t) =\Ti^ z^(t) - \lK2Z2{t) = ^^K^z^^{t)sgn[z^^{t)]- jK2Z^.2(t)sgn[z^2(^)] 

(3r53) 


The correlation function of z is 

E Iz(t^)z{t2) }= K^E (z^(t^)z^{t2)} + K 2 E (z2(t^ )z2(t2)> 

- E {z^{t^)z2(t2)} - E {Z2(t^)z,(t2)} 


{3-54) 


These terms will presently be computed 


El = E {z^(t^)z^(t2)} = E {Zg^(t^)z^^(t2)sgn[z^^(t^)z^^(t2)]> (3-55) 

Conditioned on d and d this can be factored 
11 21 

^l"^d^^ ,d2^^^^^sl ^^1 ^^sl ^^2^ 1^1 1 ECsgn(z^^ (t^ )z^^ (t2)) ,d2^]} 

(3-5£) 

Using the previously defined function one can write for the first 
conditional expectation; 

E[Zsi(t^)z^l(t2) I d^pd2^]=[T^P2d^^+a^]A^(t^,t2) (3-57) 

for the second expectation one can write for A-j=l 

E <ssn[z^,(t,)z_.,(t2)] }= E {sgn[z^,(t,)]) = 1 (3-58)' 

hence one finds 

^^1^2^ <^21^^ ^ A^(t^,t2) (3-59) 

The same method can be applied to the second term ^n (3-5^. The third 
term can be factored; 
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E3 ^ E {z^(t^)z2(t2)} = E {Zg-j{t^)z^ 2 (t 2 J^ 9 n[Zcl(t^)Zg 2 (^ 2 ) 3 > 


= E {z^-|(t^)sgn[z^ 2 ^^z 33 ^ ^ ^ 

Let E^^ be the first factor in (3-60) 

E31 ^ E {Z2^{t^)sgn[Zg2^^2^^ ^ 


= -I, ^[P^ E{d2i(t^) sgn[-T2 d22 (t2) + ^52(^2^ 3 


+ E{ Ng^(t^ ) sgnC-Tg ”4 f^ «^22 ^^2^ ^s2^^2 ^3 

t r-^2 ^ ^^22 ^^2^ 

' -Tl <1,2 j<“2l(‘l) -=>-f 


+E„ H 

Ns^,d22] si 1 


^^2 

r- TgvI^g ^22 ^^2^ + E{Ng2l Ng^ > 


where 


E{N52|Ns,)= E{Nj2t*2> l'*Sl(‘l> ’ 

Var{Nj2 I Nj,) = E{ [N52(t2> ' - E < N52 1 N^, ) ]^ I Nj,(t,) ) 

I 

For convenience we define the function 

| 1 ; kTg-Cn-l )T-j < t^ < kl2+T^ ;kl2 < tg < (k +1 )'^ 2 * 

B(t^,t2) “j 

(O; otherwise 


(3-60) 


(3-61) 

(3-62)' 

(3-63) 

k 

(3:64)- 



Then for 


JZinCL 


om 


li(t^,t2) = 0: ^sl ^ ^ °2 


(3-65) 


SUi.y = 1 : E( "s2 I "si > = " 3 , . Var(N32|N3,)=Si„E 


since N ^2 of n independent samples of Ng](t)* Hence one finds 


-31 


={t, F, 


erf 




JTo. 


2 a, 


^^ 2 ^f^ 


exp 




to) (3-66) 


The second factor in (3-60) can be computed the same way and E 3 is found to be 




Ej erf 


T2'T2 1 


>]T a 


2 J 


2 a^ 

— - exp 
\12 IT 


■ 6)1 




L'TF.,, 


2^2 
N2 irn 


6 )!«. 


.U) 


(3-67) 


Finally, the last term in (3-54)' can be found from by interchanging subscripts on t. 
Substituting E^ through E 4 into (3-64), averaging over t^ and integrating over 
(tg-t-j) gives the noise spectral density at zero: 


15 




dgk"^ 
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ont 


^0 2 
ir " 


(i+^“ 


1 

I 


+ n X 


(- 


Vi 


1 

+a) 


1 

I 


) 


-2 


?1 X 

TW 




exp 



^1 

iray 






(3-68) 


where we used the optimum gain ratio (3-13). The phase noise 

variance and squaring loss can now be found from -(3;;^31) and (3-37), respec 
tively. 

3.5 Performance Comparison 

In this section we will compare the performance of the 1m — 
pi ementations investigated in the last section. In addition we will 
investigate how the power unbalance affects the performance of the 
carrier recovery loop. In Figs, 3.4 and 3.5-. squaring loss is plotted as a 
function of the power split for balanced and unbalanced (10:1) data 

rates, respectively. Clearly, the implementation with hard-1 imiters 
performs better for most conditions and is the only one 

which is able tO' track the carrier over all data rate and power 
splits. To investigate the effect of power unbalance on synchroni- 
zation we assume that the bit energy in both channels is the same 

P^Ti = P 2 T 2 . (3-69) 
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Assuming a two channel Costas loop with hard-1 imlters^ the squaring loss 
is rniim! from (3-37) 



[erfjc-^,1^ exp(-c)]^ 

5(l+x) + l-2>|ABrerf(xrr)+ ^ exp{-5)]^ 


1 + 2 ^^ + ^ 
20 +^^) 


(3-70) 


where approach theBPSK 

case for x-»- 0 since Ejj/Nq is constant. To find the BPSK performance n is 
kept constant in (3-68) and \ is set to zero). Eq*(3-70) J^s plotted in 
Fig-3'.6as a function of A. for different values of C. Clearly, balanced QPSK 
is optimum for most values of C. Only for low values oft the optimum is moving 
to low values of A. 



1 


in K^oyn 



Figure 3V6. Squaring Loss for Different SNR, 
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3.6 Unbalanced QPSK Carrier Reconstruction Loops ’Motivated by 

MAP Estimation Theory 

Fig. 3. 7 illustrates the UQPSK carrier tracking loop which is 
motivated by the MAP carrier reconstruction loop as derived previously 
Here the integrate and dump circuits have been replaced by low pass 
filters (LPFs), the digital filters in the two arms are combined 
into a single loop filter and the variable phase oscillator is 
replaced by the VCO to accomodate channel frequency uncertainty. 

The nonlinear elements are G{x) = tanh(x). When the signal-to-noise 
ratio is small then the nonlinearity G{x) can be approximated by 
G(x) = X and the loop is recognized as two Costas loops in parallel. 

As any system in which the same signal is processed via different 
filters the problems of arm imbalance and mismatch filters must 
be contended with in the hardware. On the other hand, at high signal- 
to-noise ratios the hyperbolic tangent in Fig.-l can be replaced by 
the hard limiter G(x) = sgn x. This carrier tracking loop for 
UQPSK is essentially two polarity (modulation wipeoff) type loops 
in parallel. It functions to wipe the modulation off the inphase 
and quadrature channels thereby creating an error signal for 
controlling the VCO. In passing it should be pointed out that 
there exists the modulation/remodulation type of carrier tracking 
loop which is similar to this mechanization; however, it is believed 
that its hardware implementation is more complex and it will not 
be indicated here. To match the loop implementations as closely 
as possible with the solution of the MAP approach, we choose the 
gains C-jjCg such that they compensate for the DC gain of the integrate- 
and-dump filters (assuming unit DC gain for the passive arm filters) ,i.e. , 
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Figure ^.7 UQPSK Carrier Tracking Loop Suggested by MAP Approach. 














oyn 


c. = K.T. i=l,2 (3-71) 

where T^. is the bit time in channel i and K. are the gains obtained 
from the MAP approach. K 2 /K^ was found in [3-4] 

<y^i>opt = (3-72) 

The optimum gain ratio for passive arm filters is therefore 

(Vf^l’opt ' • (3-73) 

3.7 Performance Analysis of UQPSK Carrier Reconstruction Loops with 
Unknown Bit Timing 

3.7.1 Performance Analysis of a Two Channel Costas Loop with 
Passive Arm Filters 

This implementation is illustrated in Fig.3.7 with G(x) = x. 

The model for the received signal is defined by eqs.'(3-l), (3-Z) and 
(3-14) of Sec, 3.2. Neglecting double frequency terms the correlator 
outputs of Fig.3'7 . 


6^ (t) = dj(t)co3 cp-h qj + N^(t) cos cp 

- N^(t) sinqj .-(.3-74) 

~ dj(t) sin cp - ^ 

+ N^(t) cos 9 ' (3-75) 
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om 


where cpss 0-8 and 9 is the phase of the VCO output. Defining 
the filtered signals and noises (using the Heaviside operator) 


= G.(p)[d.(t)} 

(3-76) 

N (t) = G (p)tN (t)l 
Cl 1C 

- (3-77) 

N .(t) = G.(p){N (t)} 

Ol 1 

- . (3-78) 


the signals at the arm filter outputs can be written as (leaving off the 
variable t for concise notation) 


^ ^2i ^ci ^ “ ^si 

"(3-79) 

z , . = G. tyP7 'ii • 'P - . cos q? + N . sin <p + N . cos 

li 1, 1 ii ^ Z Zl ^ Cl SI 

(3-80) 

For the multiplier outputs one finds 




-(3-81) 


where we introduced the noise 


N.(t.E,) ^ dj.N^. t VP', d,.N^. . 


2 isin 2p 


cos Effl+ d, .N . + */l^d N .] sin 2p+ [n .-N ^ 

^ 4 Cl ^ 2 2i si-* Cl SI Z 


' (3-82) 
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Let now 




sin Zcp 




(3-83) 


+ Nj(t,2<p) . N2(t',2<p) 


Assuming zero VCO instabilities and VCO timing voltage, the 

A 

instantaneous phase 6 of the VCO output is related to its input 
through 


e = 


K^F(p) 


(3-84) 


Hence for a constant signal phase 9, one finds 


pcp= -K^5'(P)([C^{Pid^rP2'^2l^-'^2<^l'^U-V22)l^ + 

I ~ (3-85) 


where 


cos 2(p sin Zcp 

N^(t, jp)=:2[N^(t, cp)-N2(t, cp)] 2 ' ^ W 2 


and 


(3-86) 


X(j - 2[c^yPp^ ^11^^21 ’ ^2*^^1^2 ^12^22*^ (3-87) 


x_ ~ C^(N^ -N^ ) - C^(N^ -N^ ) 
^ si’ 2^ c2 s2^ 


(3-88) 




(3-89) 


" ^^^l^'/^i^U^cl‘*‘*^‘^21^sl^“ ^2^"^‘^12^c2‘^*^2'^22^s2^^ 
A . 2 — ^ 2 _ (3-90) 

(3-91)- 
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All these variables have mean zero. Using the independence of 
data and noise it is easy to show that they are uncorrelated, hence 


the correlation function of N is 

e 


sm 2(0 


cp)N^(t+T, qs)|<p} = + ^ 






Asstiming the correlations 


i = 1,2 

si ci 1 


(3-93) 

(3-94) 


it follows that 


Rj,2(t) = (0) + 2R^ (T) i = 1.2 

R„2j,2(t) ^ R„ (O)R^, (0) + 2rL(t) 


V'Z 


N, ' ' N. 


NN' 


(3-95) 

(3-96) 


and the correlations of the x-processes are 


-cfc^R, , (t)R (t)] 

^ ^ ^ 11^12 ^ 21^22 

Efe^x J = E J = 4(C^[PjR (T) + P^R^^ (t)]R„ (t) 


(3-97) 


11 


‘El 


+ C^CPiR,^^(T) + P2Rd^^(T)]Rj,^(T)-2C^C^Cp,R^^^^^^(T) 


+ A 

^ 21 ^ 22 ' 


(T)]Rj^(T)} (3-98) 
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(.3-99) 


Substituting (3-97) ,(3-98) into (3'-92)rs.ettingcos 2cp= 1 and integrating 
over T yields the noise spectral density at u)= 0: 


N ® A * 

-=2., p P {c^f R, (t)R, (TldT+cM R, (T)R^ 

2. 1 21 IJ^^ d^j dj2 *^22 

(T)dT 



•fo 11 bX 


+ J* !>*l«d, <"> ^2 V *"n2<’-'1 '^N2<^>'*’- 

-«o 12 22 


■ - ^<=^2 f t^l^d. ,d, (■^1 + P2«d„d„f^>+W>1^NN<‘'>'*4 

11 1a a1 

(3-100) 

We note that 


® - .o__- 

«m<^> = I -# 1S'«I 

-<a 

(3-101) 


(3-102) 

(t) = J S^ {f)|Gj^(f)l^ e^^^'^df 

• ik -«» i • 

' (3-103) 

^d,d <^> = r 

il i2 -« 1 

(3-104) 


Substituting these relations intc (3-100) yields 
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om 


n; Ng rN 


r [ X * ^2®, 

^ 1 


2P P 

(1) t- jj 'S (f)S (£) 
"o d 

£t 


■] 


.[CjlGjWf - cl\a^{f)\^f if 
We linearize the loop equation 


(3-105) 


•pcp= -K^F(p){Cc,(P^d^j-P2d^,)-C2(P^o/^2-P^d2^)] <p+N^(t,2cp) 

- - (3-106) 

and replace the qxiantity in brackets by its mean value. Then one 
finds for the phase error as a fxmction of the noise 


KyF(p) 

'■^■■p+cfCPiR, (Ol-PjR, (0 )]-c^CPiR, (pj-p^R^^^to)] 

11 Z 1 1 m 


N (t, ■«?) 


‘22 


or, setting 


(3-107) 


s»( 0 )'=* C,[p,R, (O)-P-R, (0)]-C^[P R (0)-P,R, (0)1(3-108) 

^ ^ ‘^ll ^ ^21 ^ ^ ^12 ^ "^22 


The phase error variance can be found from [3-4] 


2 


a 


«P 


[s>(0)f 



(.2-109) 


where S’(0) is the slope of the loop S-curve at «p= 0 and B_ is the 

L> 

loop bandwidth. The performance of a single Costas loop tracking 
one of the channels can be found by setting or to zero. Using 
the optimum gain ratio (3-73) and se-tting 

■ ^2 = = il: 1 ^2 = ^ 

One can write for the noise spectral density 
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N’ N.P . « 

0 0^ ^4 p 

2 ’2(1+X) ■ J 


Nq(1+X) 

”'2P 


+ S, (f) + >3^ (f) + 


2\P 


S, (f)S^ (f) 




[ |Gj(f)l^-nX|G2(f)f-]^df 


(3-m) 


and for the a lope of the S- curve 


S«(0) , = 


Cjp 

x+x 


{R, (0) - XR, (0) -nX[R, (O)-XR, (0)]} 


(3-112) 


‘11 


‘21 


‘12 


*22 


Ab in the case of active arm £ilteraii3-4] thi? implementation does not 
txacJc on halanced QPSK since for n s 1 the adope ia 

0)C1 

11 #, 

which is zero for X= X. 

. s 

3;7.2 Performance Analysis of a Two Channel Costas Lioop with Passive 
Arm Filters and Hard- Limiters 

Since this loop is identical with, the one analyzed 

in 6.1 » except for the hard-limiters the demodulator outputs must 
still be characterized by (3-l5K3-16).The filters are also unchanged 
from the last section; however, to keep the problem manageable the 

/ 

effect of the filter on the signal will be considered in this section as a 
reduction in signal power only, i.e. the filter outputs are still 
square-waves corrupted by Gaussian noise: 


- C. cos qj + ^/P^* ^2 cp + N . cos <p - sin cp] 

^ ^ ^ cx a (3-114) 


dj^ sin qj - *^2i ^ ^ci ^ai ^ /(3-115) 


Using the same approach as for the active arm-filter case [3-4] we 
find the slope of the loop S-curve at <p= 0 and the noise spectral density 
for qj = 0, 


page is 
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C 


3- 7=^. 2.1 Loop S-Curve 


where 


Setting 


is defined by 


S(<p) = E(el<p} 

(3-116) 

C = Cj - 62 

(3-117) 

6j = 6^1 6gn(6^l) 

(3-118) 

^ = V 

(3-119) 

s.fe) = teiltp} 

(3-120). 


we find 


Sj(cp) = E{z^j sgn(z^^)jcp} 

= sin cp- */^21 ^2 ^ci ^sl 'PI* 

‘ agnCv^,d, cos cp + */P^ ^2 'P ^cl ^ “ ^sl 'pII^^ 

(3-121) 

Conditioned on dj,d^ this escpectation can he factored since the sunis 

of the two noise components are again zero mean independent 

2 

Gaussian noise processes with variance a . = N„B. where B. is the arm 
. 1 0 1 1 

filter bandwidth 

Sj(cp) = ^ sin cp - .Jp^ d^ cos cp + sin cp 

2 (3-122) 

+ cos cpjd^, d 2 ,cp]E[sgn(VP^^‘ d^ cos cp + */P^ d^ sin cp 

+ cos cp - sin cp) |d^, d^, cp] !cp} 
the conditional expectations are easily computed; OK[CjIN‘Al3 PAGE IS 

DJI EOOR quality 
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cm 


Sj (cp) = I |;yp^ d^ sin cp - d^ cos qj ] 

[ J^ll <^1 *=03 'P + ^2 



THs is readily averaged over the four pairs (dj^^d^) 

( f ^ *’■ *^^21 

Sj(cp) - iCj|(^sincp cos €p)erf|^ " ~ J2Y ^ J 


■]l<pj 


(3-123) 


[ a/P,/ cos o - a/P,! ^1) 

I 

(3-124) 


Similarly one finds for S 2 (cp) 


SjW' = iCj |(yPYJ cos <p sin <p)erf ■ 


a/Pj 2 “VP^ cos cp 


22 


1 


+ Cy^^COS Cp - VP22 


a/^02 

VP^ sin cp ^ cos cp Tj 


(3-125) 


A/Z'a^ 

To find the slope at cp * 0 we take the derivative of and with 

* i 

respect to cp : 

5« 

■ A^- Ai j_V* c^j i.-*- V" 2a** 

(3-126) 


I 

- f /S-- 2 I 

3^(0) = C^|a/^ " *^1 a/I?®"^ } 






'"*^22 


+ 


The slope of the S-curve at <p = 0 is then 

\ 

s*(0) = S^(P) - S^(0) 


2 ^22^ I 

^ (3-127) 


(3-128) 
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3.7. 2,2 Noise Spectral Density 

The correlation function of the noise can be split into auto- and 
cross -correlation functions of the multiplier outputs e., i s 1, 2 
= E& (t)e (t+T)] = RjjjCt) + •- 

(3-lb) 

where 

(t) = E|:c.(t)e.(t+T)} (3-130) 

X X 

(3-131) 

The four terms of (3-129) will presently be computed. For the first 
one finds 

. Cj^E[(-^2i “^2 + (3-132) 

.ECsgn[(7F^ dj + N^i)C./F7J dj^+ 

r 

j 

since iiidependent. The first e:q>ectation is simply 

E{(-.^d2+N^j)(-7FJJ'd2^+N^^^}eP2^R^ (T) +Rj^3i(T) 

2 

The second expectation is (3-133) 

E{sgn[(#^ dj +N^i)(.A^dj^ + (3-.134) 

= dj Idj- di,l) 

Calling the conditional expected value M(dj, it can be shown 3] 
that 
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JZmd 


JP ' 


« P?(T) 




(3 


0 j n=l 


where He is the Hermite polynomial and Pj^(t) is the normalized 
correlation function of the noise 


Pl(T) = — 


(3 


Averaging over d^jd^^^then yields 


^-1 d = 3[R^ (t)+1]m( 1,1) +i[l-R, (t)]M(1,-1> 
“i°1t 1 1 


and becomes 




(3 


*11 \ ■„ 2 


0j n=i ^ ^ 

In the same way one finds R -(t), Finally, one finds for the 


cross-corr elation 




om 


- 135 ) 



- 136 ) 


- 137 ) 
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= CjC2E[(-,yp^d2 + N^j)agn(-VP^d2^+ 

.E{(yP^dj^+N^^^)sgB(^ dj +N^j)} 


(3-138) 


Conditioning the first expectation on d_, d, , N , it can be factored 

u 5i 

Et (-^ d^ )sgn(-VJg ■i 2 T+N 32 ^ Idj, N j} 






’^sl v^'vKTTi?) 


•where 


m(N.,) = ECN^^^lN^jl = 


ri2(T)o-2 




IZ <7ja2 


(3-139)' 


:r"^[S^(f)GJ(f)G2(f)] ■ (3-141) 

X 2 


(where ^ stands for the inverse Fourier transform) 




(3-142) 


Hence the first expected -value in (3-136) is 


E ^2 + 1 t"'' 2 

Interchanging subscripts one gets the second expected value in (3-138) 
Substituting correlation function 

of the noise 
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(t) 


“ A i 

2 /■ 2 V 

• - exjjf- He , f ii 

^ ^ a^^n! n-lU^ . 

- 


I 1 

~)i&+K^ (t)) 


- (-l)“(l-R^ (T))]} 

' 2 

+ ‘^zl^U^dj <■'> + ®Nc2<''>^ ’ 

^2' / **22N^ '2 /-^N 


-Ci‘=z{vPiTV(f)«f[;jXj Ojrjj(i)«p .(— ^)} 

2ff2 

+yl'’2’'l2<-’-)==^ -(^)) 

* «- 2 J 2 q 2 

r* '*^ * 

+7f a^rjjMe^cp -(— y^j. 
*- --i-i 2o, 

I 


-.(-i)”(i-a^ (T))]} 
3 

■j2(T)e==p -(-^)} 

2c, 


wliere 


(3-144) 


(T) = 


( ^ I I ^ 

0 |tI>t. 


(3-145) 


“ ^ci'’'> = ^’^[-r l°iWl^] - (3-146) 


rj 2 (T) = 


1 ,.:ir”^o * T 

— 3 [ ^ GJ(f)Gj(f)J 
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The integral of over t is the noise spectral density at f = 0. 

To simplify the notation we define the arm filter bandwidths B. = aR^ 

and introduce the parameters = P/(l+X), = XP/ 

P, . = P.S, , (G-/C, = nX, p, = PT, /[2N a(l+X )] and the notation 

11 i d.. 2 1 i 0 

I{x) = j’_^ x(T)dT (3-148) 

Then the noise spectral density becomes 


2 f r A 

' (l+X) ) H V3n) 


^n)"^ T^P 


" ■ 

erl (*/psl ') 

J 11 


+ Xb 




(1+X)I(H^^,2^^2^ 

T^P 


erf (;^Xnps ') 


‘21 2 


“ exp(-2ps ) 23 


^,„2k „ „ , 

^^^Nsl^dl^d2^ 


"11 k=i L aj (2k)t 

, ^<tl'^d2) 2 

1“ _~2k _ 




nAs r 

®12 2 _ 

+ -V; ^ exp(-2pnXs ) 

1 °22 k=l 


^^^Nc2^dl^d2^ 


" ^22 k=l L cr^ (2k)! 

■ r 2k+l 

• I{r2^ ) 1 
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txj (2k-l)! ^11 
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- Zxyn'fl erf(Vps, ')erf{*/piiX3 ')V^ sT ‘ 

^ ^ll ^22 ^12 *^21 



*22 


‘11 


erf(Vp3 , ‘)exp(-pn\s )I(r.^R,,) 
®il "^22 




(3-149) 


Using the same notation in (3-128) the expression for the slope of the S-curve 
takes the form 


+ erf(V5S^) - 2 /^ e:=p(.pnxs^^^)3^^^} 

and the squaring loss is 

No[S>(0)f 

*L " N' P. {3il51) 

0 - 


3.8 Performance Comparison 

The squaring loss has been evaluated and plotted in Figs. 3.8 and 3.9 
for balanced and unbalanced data rates for the following loops; 

A. Two channel Costas loop with active arm filters (analyzed in [3-4]). 

B. Costas or squaring loop tracking either I (Bl) or Q (B2) channel 
signal works either for X « 1 or X ^ 1 (depending on which channel is 
being tracked) but not for both 
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C, Two channel Costas loop. This loop coihbines the power of I and 

Q channel and works therefore over all power splits for unequal data rat 
None of the loops above can track a balanced QPSK signal. 

D. Two channel Costas loop with active arm filters and hard- 
limiters (analyzed in [S4] . This implementation performs well 
for all data rates and power splits but has several false lock 
points [3-4]. 

•E, Two channel Costas loop with first order Butterworth arm 
filter and hard-limiters. This loop shows an infinite phase 
error variance for X«. 1 and unequal data rates. This is due to 
the large filter bandwidth chosen (3 times the data rate) to satisfy 
the assumption of undistorted pulse shapes. This departure from 
the design suggested by the MAP approach invalidates the choice 
of the channel gain settings and performance can therefore be 
improved by numerical optimization of the gains. 

The fourth power loop is another implementation capable of tracking 
a balanced QPSK signal. It is nevertheless restricted in its use by 
self -noise which prevents carrier tracking for X , 2 [3-7], 
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4.0 SHUTTLE LINK DEMODULATOR ANALYSIS OF SQUARING CIRCUIT MECHANIZATIONS 
ITS ED IN SQUARING LOOPS 


,4.1 Introduction 

It is well known that the Costas loop of Figj4- la and the squaring loop 
of Fig- 4-Tb have the same theoretical noise immunity in both the acquisition 
and tracking modes. However, in the implementation of a squaring loop, 
mechanization of the times two multiplier is an important consideration inso- 
far as system performance at low signal -to-noise ratio, is concerned. Con- 
siderations which must be accounted for in the choice of the "squaring" 
approach include a wide dynamic range with respect to the signal level, good 
thermal stability, and accurate square law response over the dynamic range of 
input signal and temperature levels of interest. The performance with an 
analog multiplier usually indicates degraded signal -to-noise performance 
relative to theoretical. In an attempt to overcome the degrading effects 
at low signal -to-noise ratios, an alternate approach to the implementation of 
a squaring circuit was considered. The alternate mechanization of the loop 
is illustrated in Fig. 4-2. As far as the authors can- tell the tracking 
performance of this limiter/multiplier implementation of a squaring circuit 
incorporated in a PLL has not been addressed in the open literature. 

The purpose of thissection is to investigate the performance of this 
implementation as compared to the known performance of the ideal squaring loop 
case given in [4-3]. In particular, relative comparisons are made on weak 
signal suppression factor on the phase error signal in the second imple- 
mentation, error signal SNR as a function of the input SNR, effective loop 
SNR and squaring losses as functions of input filter ban^dwidths and input 
SNR's. From this computation the degradations on loop performance due to 
this limiter/multiplier implementation from the ideal squaring loop case 
is established. 
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Figure 4 1 








4.2 Preliinlnarles 

Fig. -4»-2 shows the limiter/multiplier type of squaring loop under 
consideration. In that figure, h-j(t) is a bandpass filter with center 
frequency and equivalent single-sided noise bandwidth B^(B^ « uq), 
given by 


B 


i 



T )l ^ 

*0 |H(>q)1^ 


dm 


(4-Jj 


where H is the transfer function of the filter. The input signal to the 
limiter/multiplier combination is defined to be of the form 


with 


x(t) = s(t) + n(t) 


s(t) =41” A sin { upt + e(t)) 


(4-2) 


(4-3)- 


where e(t) is the information bearing signal, and n(t) is the narrow band 
noise represented by 


n(t) = [n^(t) cos oj^t - n^{t)sin Upt] 


(4-4'): 


where n (t), n (t) are zero mean uncorrelated Gaussian processes with 

V* ff o 


E[n^(t)] = E[n^(t)] = 
=■ E[n^(t)] = NjB, 

E[n^(t)ng(t +t)] = 0 


(4-5-)' 


In (4-5)N is* the one sided spectral density at the input to the BPF 


om 
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The noise process n{t) xan be further written as 

n(t) = N^(t) cos (cj^t + 9^{t}) (4-6) ) - 

where N^, 8^ are the noise envelope and phase respectively given by 

N^(t) = ^ n^U) + n|(7) 
and 



Combining (4-3) ,.('4-6j. through some trigonometric identities we can re- 
write the input signal x(t) in the form: 


x(t) = /S' [(A-Mg(t)) sin 4'{t) + N^{t) cos *(t)3 (4-7) 


where 

♦(t) = u^t + e(t) 

N^.(t) = N.(t) cos (0^(t) - e(t)) (4-8) 

NjCt) =■ N-(t) sin (9^(t) 8(t)) 


The noise processes involve the information bearing signal 9(t). 

To obtain the statistical properties ofN^,N^we can rewrite (4-8) as follows; 


NjCt) 


n (t) cos 0(t) + n^(t) s1n e(t) 
n (t) cos s(t) - 0 (t) sin e(t) 


(4-9) 
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which gives directly the following, assuming the noise process and the 
information bearing signal are statistically independent 


E[N^.(t)] = E[Nj(t)] = 0 

ECN^]= E[N^] - a^/Z (^-10) 

E[N^(t)N5(t)3 = 0 

If in addition to (4~$ we.- define the normalized input noise autocorrelation 
function to be 


( t) = r {T)coswnT . 


'1 


n n 


-e><-f<a 


where r^( is low pass and has the properties 


'•„(0) = 1; |r„(T)|<l; f dr = ^ 


(4-11) 


then the auto- and cross- correlation functions for N^{t) and Nj(t) can be 

2 

written in terms of r (t) and a as: 

n n 


ECN,N^,; EfN^N„] = r„(r) cosie^ 

eCN.N^^] - -ECN^N^^] r„(r)rnse7 


(4-12) 


where 


ae^=e(t)- dlt + r) 


^lOd- 



Xi„C 


om 


4.3 Limiter Output Signal 

The ideal limiter function is defined such that 


y(t) - sgn(x(tl)>r' If x(t) > 0 

'-1 ifx(t) <0 

The input process x(t) can be written, from (4-7), ^as follows 
x(t) = /~2 v(t) cos (-}{t) - -y{t)) 


\ C4-13) 


(.4-14) 


where 


= i/{a-n/+ N 


Y s tan 


-1 




(4-15) 


Using the complex Fourier transform approach of relating the limiter output 
time function to that of the input, one obtains 14-]] 


y(t) » 


rs(j 

2irj J 


G(jx) exp [j\ x(t)] dx 


(4-16) 


1 _ f ^ 


exp [jx v(t) cos («(t) - y{t)) ] 


with G(jx) = 2A t the transfer characteristic of the hard limiter [4-T],^-and 
C the dynamic path of integration. Expanding the exponential function in (4-16) 
by the Jacobi -Anger formula 


jzcos 0 - ^ 1^(2] cos ke 


k=0 


^ StMEEQJ 


-no- 
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where 2 for all k > 1 and = 1 » and where are modified Bessel 

functions, the limiter output y(t) can be integrateu to be [4-2]: 


y(t) 


^ cos {(2k <• 1) C «(t) - Y(f)]{ (4-17) 


It is clear from (4-17). tHa-t’. only odd harmonics are present in the limiter output. 
^♦4 Squarer Output 

To obtain suppressed carrier tracking, the data modulation Sit) 

[bi-phase modulation assumed] has to be removed to create a sinusoidal signal which 
is tracked by the PLL. This is accomplished by multiplying the limiter 
output (a stream of - I's) by the incoming signal, and then obtaining the 
hamonic' around 2tu^ by passing the multiplier output throught the zonal 
filter h 2 * The multiplier output can be written as: 


x(t) * y(t) = /Tv(t) cos ($(t) - Y(t)) 


(4-18) 


4 * 

I E 

k=0 


cos { (2k + 1) [$(t) - Y(t)] } 


The second harmonic term can be selected from (4-18) to be the following: 

z(t) = v(t) cos [2 ($(t) - Y(t}) 1 (4-19) 

Note then that in the above discussion the filter h 2 is assumed a mathematical 
entity which selects only the second zone component. In reality, a physical 
bandpass filter can only approximate this condition since the spectrum of 
cos [2 (4>(t)-Y(t)] extends on all w; however, because of the assumed narrow band 
approximation (i.e., B-|«wq), the error is small. Notice also from (4-19) that 
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since bi-phase modulation is assumed (i.e., e(t)= 0 or u), z(t) is actually 
a sinusoidal signal at twice the earner frequency Wg when the noise effect 
y(t) is negligible. 


4.5 Weak Signal Suppression 

To obtain the loop error signal, the squarer output z(t) is mixed 
with the local reference signal 

r(t;^) = -■/2' sin 2[w t + <}»] (4-20) 

» 0 


where ^ is the phase difference between z{t) and r(t;ij)). Onitting terms 
the error signal that the PLL tracks is then 


6(t) 




- (4-21 ) 



+ 


VNj + (A - N J 


2 N . (A-N.) 

sin 2a •+■ cos 2<i» 


It is of interest to compare this error signal with that of the perfect 
squaring loop (See [4-3], pg., 59, Equation 2-69.) 


^sqt^) “ sin 2 ^ + [2 N^(A-N^)] cos 2iJ' | (4r22) 


in which Kp are rms VCO output in volts and multiplier gain respectively. 

Disregarding the multiplicative constants both e(t) and e (t) can be written 
in the form ^ 
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0(t) = a sin 241 + 6 cos 24 


e (t) = a' sin 24 + 3 ' cos 24 


(4-23) 


where a, s, and a‘, 3 * are random processes as defined in .eqs.(4-21) and (4-22) 
respectively. For the perfect squaring loop case (eq. '4- 22) the 
following equations are readily verified from (4-10) 


E(a') = E[(A-N^)^ - N^] = 

E(S') = 2.E[N^(A-N^)] = 0 


'(4-24) 


Thus for the perfect squaring loop the mean error signal e (tj is 

2 

directly proportional to A sin 24 , and the terms' cos 24 can be 
considered as a zero mean equivalent noise process. Notice then in 

I 

the perfect squaring loop case the mean amplitude of the error signal 
e(t) depends only on the input signal amplitude A. In addition, the 
SNR of the loop error signal can readily be computed from the first 
two moments of e 5 q(t) to be the following 


SNR 


Ejq Var(3 ' ) 


~2“2~ T 
2A 

n n 


C4-25) 


This can then be written in terms of the input SNR as 

2 


SNR 


'sq 


1 + 2p. 


(4-26) 


where p. is the input signal to noise power ratio, defined by 


N B 
0 1 


(4-27) 
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The output to input SNR rations given by 


SNR 

Sq _ 

" V + '2p^‘ -■ (4-28) 

which approachs 1/2 as and approaches p- as p.— ►O. Fig; 4 t8 

shows a plot of this ratio. 

However, in the case of the limiter-multiplier implementation of the squaring 
loop, there is signal suppression on e{t), which is a fCinction of input SNR p^. . 
Suppression in e(t) will affect the loop bandwidth and tracking performance. 
However, if this can be compensated for by some control of the multiplier or 
loop gain, then signal suppression on e(t), by itself, will not affect loop 
performance. The actual performance, of course, will depend on how much the 
noise power is suppressed relative to the signal power, which is discussed 
in the next section. To evaluate the SNR^ for e{t) it is necessary to compute 
the first two moments of the terms a, 6 defined in eqs. (4-21) and (^4-23). • 
as the first two moments of a, 3 are concerned N and N occur at the same 
time and thus they are independent Gaussian with variances /2, and the 
evaluations of the expectations are reasonably straightforward. The results are 


E(a) = E 


(A-Nj)^ - n/ 




+ (A-N^)2 


= |-/TA/p. e *^i P[|-;3;p^] 


E(e) = 2.E 


" N^(A-N3) 

+ (A-N, 


C4-29) 


' C4-30) 


driqetah page b 
OM EO.Q8. fiUALOT 
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Figure 4.3. Normalized Signal to Noise Ratio of Loop Error Signal 
as a Function of Input SNR. 
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E(a^) 


= E 


(A-N3)^ - 


ijnl * (A-Nj‘ 


= 5^!- 0 + Pj) + g-p, e'^l Ft4; 6; p,] 


(4-31) 


E(r) > 4-E 


“c 


V“c ^ 


2 2/ P 2 ) 

= ^-|2(l4.p.) F C4; 5; p^]|(H2): 


where F is the confluent hypergeometric function [4-6] defined by 


E [a. Y» x] 


y (g; 1; v) 

^ V I f V • 1 • 


v=0 


v! Cy; 1j v) 


w1 th 

(m; d; v) = m (m + d) ••• (m + (v-l)d) 


{4-.33)., 


and 

(m; d; 0) =» 1 


For §“33) to converge-, , in general, it is required that jxj <1 and y 5* 0> 

-1, -2, .... . However, for the cases of interest in Equations (4-29) to -(4-32) 
the particular confluent hypergeometric functions involved can be written 
in terms of simple functions and modified Bessel functions (See [4-4] ,p. 1073) 


F C f; 3; 3 



(4-34) 


F C 4; 5; p- ] 


(pj - 3pf + 6 p. -6) + s] 
•"i 


(4--35) 
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Written in such forms then it is easy to see that the convergence of (4t29)-(4-32) 


actually does not depend on the value of p^. 


Since for large p^/2 


Pl/2 


t,n 


thus the mean amplitude of the error signal c(t), as is given from 

(4-21) (wtth the multiplicative constant 4/3 neglected). 


lira lira. 

E [e(t)l = 

0 .— *<» ^ ^ Q .— ♦« 


(A-N^)^ - n/ 




sin 2iji 




11m 

p, - 


|/7A/f: fer" 


P4 r 


Pi I P4 


sin 2({> 


(4-36) 


A si n 2(f) 


The weak signal suppression factor on E [^{t)] as a function of p^. can be 
given as 


g- [£( t)] . 


^ j — „ 
’ r Q ^ e 


Pi 

T 




<’ ' 


(^) ] sin 2^ (4-37) 


The multiplying factor on sin 24* in eq.._(4-37) is illustrated in 
Fig-.-4'.4, as a function of the input SNR, 
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Figure 4.4. Signal Suppression Factor on Loop Error Signal As A Function of 
Input SNR. 





Similarlly to eq. {4-48}' we can compute the ratio between the SNR 
of the error signal e(t) and the input SNR The following result 
follows directly from (4r34) and (4-35) 


SNR^ 


(f /7 e-°1 F[-4; 3; p.])^ 

5 jan to,) - F[4; 5; p,] | 



This ratio. is compared to that of the perfect squaring loop in Fig* 4.3* 

It is of interest to observe the following limiting cases for the above ratio 


SNR , 

^ ^ as p.| — ^ « [4-39) 


SNR 


9ii 

3? Pi Pi 


0 


(4-40) 


As far as this ratio is concerned the limiter-multiplier implementation is 
equivalent to the perfect squaring loop as p. -#*«and is ~0.54 dB (9V32) 
worse than the perfect squaring loop as p.j —*■ 0. 
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4.6 Effective Loop SNR and Squaring Losses 

As far as loop performance is concerned the essential factor is the 

effective loop SNR { [4-5] ) which depends on the noise power spectral 

density N_ of the equivalent noise process in the error signal c{t)» 

6m 

around zero frequency. N is defined to be [4-3] 


Q» 

T /. os gq 

Where is equivalent to zero mean noise processes in £{t) (eqs.‘(4~21) 
or (4-22)'). . For the perfect squaring loop case the auto-correlation 
function (t) of the equivalent noise process in (22) has been 
computed by^^ Lindsey and Simon [4-3] "to-be 

R„ (t) =» 4 [A^ (t) + ^(t)] (4-42) 

sq c c 

where R^^ (t) is as given in (4-12) . From this the equivalent noise power 
spectral ^density around zero frequency N is computed to be 

"sq ' * *^”0 ^C' 

where is the input noise spectral density and is defined to be the 
circuit squaring loss 

00 

- ■> (t)dt - (4-44) 

’0 - CO ^ 
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The effective loop SNR f^ff for the perfect squaring loop is given [4-3]-i'n 
terms of^^ and the equivalent signal-to-noise ratio in the loop band- 
width (B, ) of a second order PLL = A^/N_B, by 

U» 0 u 


f^ff (perfect squaring) 



r (4-45) 


To compare the two implementations it is necessary to compute the for 
the equivalent noise processin (4r21) and to define the effective squaring 
loss in the second implementation through eq. (4-45) by computing the 
effective loop SNR f^ff- To obtain in this case, requires computation 
of the auto-correlation function of equivalent noise term in (4-21) 



f n^-(a-N3) \ 


COS 


Assuming <j) ^ 0 (for tracking) the auto-correlation R (t) is found from 

eq 


R„ (t) = 4 E ^ ^ ^ 


eq 


^1 + (A-Nj)^ 




(4-46) 


The randcro variables N^, N^, are jointly Gaussian, zero mean, 

with covariance matrix (seeeq. (4-12)) 


_n_ 

2 


/ 1 


0 


COS 

40^ r ^(t) sin i 

0 


1 

-r„M- 

sin 


r^C^) cos 


- 

sin 

1 

0 

^^r„(T) ?TF 


r„(T) 

cos && 

T 

0 

1 / 


- 121 - 



For all practical considerations the correlation times of the input noise 

n{t) is much shorter than that of the modulation process e{t), then the 

'actual covariance matrix A of N , N , N , N is given by, for all 

c s ct s c 

practical considerations, the following: 


r 


T* 






0 

1 

0 


r„(.) 


0 

1 

0 




r^h) 

0 

1 


(4-48) 


process 


For simplicity in notation, define 

X, > , X2 = 

y, » A - "s • ^2 “ * - 


Then the joint density of x-j, y-j, X2, 


are given by 


p(x 


,.y,.X2. yj) ° 4;rr .T/2 




2|a 


1/2 


- 2 r^{x) (x^X2 - (y^ - A) (y2 - A)) ] 


(4-49) 
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Where [a] = -Jl-[1 - r^tcU- With A 0 the computation of the 

expectation in (4-46} involves quite complicated four fold integrals and 
numerical integration seems to be the only possible method of solution. 

If A = 0 (which is a good approximation to small input SNR cases), the 
the expectation (4-46) can be evaluated exactly. In terms of the noise 
envelopes and random phase angles: 


V. . + . 0, = tar.-' (^); i = 1, 2 

the expectation (4-4^. can be' computed from the following integral: 

J J dv,dV2 V, Vj e 


(4-50) 


eq 4 it2|a1 ' 


0 0 


2ir 2 it 


^ j d0^ d 02 sin 20 sin 202 
0 0 


cos(?i2-0^) 


21 Aj 


1/2 


(A-51) 


The double integral on 0^ and dg can be evaluated directly to be 


2 it 2 w I 

r r 

l°’n''n^^^''l'^2 cos(02-«ii) 

j 1 dd d<5 sin 26^ sin 2gJ exp 

! 2U|1/2 


^ 2ir2 I, 


( a ] 


With this simplification them (5-51) can be evaluated (see [4-6]) 
to be 

= 3»^.' . '’o', (' 


eq 


^ [2 - r^2(t)J^^ 


(4-52) 


X I 4 F if; f; l;x^]+ 5x F[f; f; 2; x^l 

- F[ f ; f ; 2; x^] I (4-53) 
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where we have defined 


X - C ^' 


om 




4 (c) 

n' * 


(4-54) 


and where the confluent hypergeometric function is defined through 


P [a; 3; y; x] 




v»o 


(g; 1 ;v) (s; 1 
Vi ( y ; 1 » v ) 


(4-55) 


where fo; 6; k) is definedin (4-33). The effective loop SNR can then be 
computed, with this approximation for small SNR cases, from (4-37) and (4-53) 

to be 





I (^) + 


(1 


eq L 



(4-56) 


where is the one sided loop bandwidth and is the equivalent noise 
spectral density computed from Rngq^^)* Writing as N^W./4 (4-53) . 
be written as 


(c) = g(c) 


(4-57) 


where 

g(r) 


(1 M) " , 


(2 


VT 


4FII! li U.x.^1 


(4-58) 


*S t FI^; f: 2;X^1 - 


FCf-: I ’. 2;X’]} 


and 


eq 0 


cm 

,W4 / 


g(t) dt 


(4-59) 


Then the effective loop SNR can be written in the same form as (4-45) 
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where the equivalent squaring loss j^j^can be written as 

, |»i e'" [io(r^ * (1 -^) 1,(2^)] 


J g{t)w 


,dt 


( 4 - 60 ) 


where g is defined in (4-58). 


4.7 Numerical Results - 

, As an example, consider a bandpass filter with an RC transfer 
function. The equivalent low pass spectrum for N (t) or N (t) has 
correlation functions: 


^ exp(-W,|T|) (4-61) 

Assuming signal distortion due to filtering is negligible, then the 
squaring loss for an ideal squaring loop for this Rj^ (t) is computed 
[4-3] to be ' ^ 


Jx. ■ r Y '■ (4-62) 

For the same correlation function the equivalent circuit squaring loss 
fori-the limiter/multiplier implementation can be commuted from (4-60) 
where the integration on g(t) can be performed numerically. This result 
IS plotted as a function of together with eq* (4-62) on Fig. 4.5. 

It is noted that the 'limi ter/mul tiplier implementation has relatively 
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Figure 4.5 Comparison on Squaring Losses Between Two Squaring Loop Implementations 

(BPF hi assumed to have a One Pole RC Low Pass Equivalent Transfer Function 
in Both Implementations). 





om 


more squaring loss than the ideal squaring loop for low input SNR 
cases, which is expected. However, it is interesting to note that as 0 
the difference between the two squaring losses asymtotically approaches 
«0.8 dB. 

As p. becomes large, the A «0 approximation is no longer valid. 

However, it is seen from the definition of N (t) in 
(4-21) thac 


'’eq '*■ as c« 


and thus 


FL {"C) 4R« (“C) as p. -»■ <*» 

eq c * 



as -► <a 


On 

as 


the other hand, since the signal suppression factor approaches unity 
P^ -► •, the effective loop SNR approaches, as -+• - 



A2 
0 L 



as p^-.^ ■*> 


(4-63) 


From (4-63), it fs clear that-c^ approaches unity, as was in the case of 
the ideal squaring loop case ,eq. (4-62), Therefore, we conclude that the 

loops have identical tracking performance at high signal-to-noise ratios, 
as 
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4.8 Conclusions 

Implementation of a carrier recovery circuit for a suppressed carrier binary 
phase-shift keyed signal poses the interesting question of which approach to 
take. This paper develops the tracking performance of a practical squaring 
loop in which the times two multiplier is mechanized as a 1 imiter/muUiplier 
combination; this "squaring" approach serves to produce the absolute value 
of the arriving signal as opposed to the perfect square law action which is 
required in order to render acquisition and tracking performance equivalent 
to that of a Costas loop. The absolute value type circuit appears to be the 
more practical circuit to build when such things as low signal -to-noise 
ratios, a wide dynamic range of signal level and temperature variations are 
considered. In the signal -to-noise ratio region of interest, it is shown 
that an absolute value type "square law" circuit degrades the input C/Nq by 
0.5 to 0.8 dB over that of an ideal squaring loop. This also says that the 
tracking performance of a Costas loop is better, by 0.5 to 0.8 dB, than 
that of a squaring loop implemented with a limiter/multiplier combination 
for times two multiplication. At high SNR it is shown that the tracking 
performance of the two mechanizations is identical. In addition, the beat 
note level and phase detector gain are nonlinear functions of the signal - 
to-noise ratio at the input to the limiter/multiplier. This is of concern 
as far as maintaining the design point loop bandwidth and damping as the 
signal level varies. Finally, the weak signal suppression factor is derived 
and plotted as a function of input signal-to-noise ratio. 
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5.0 DETECTION OF SHUTTLE SIGNALS TRANSMITTED THROUGH THE 

TRACKING DATA RELAY SATELLITE (TORS) 

5.1 Introduction 

This section is concerned with determining the effects of 
uplink and downlink noises' and nonlinear signal distortions on the 
detection of Shuttle signals transmitted through the TORS Satellite 
repeater exhibiting both AM to AM and AM to PM characteristics. 

The performance of the Shuttle/TDRS forward and return link is 
determined in terms of probabilities of symbol error as well as 
transition probabilities when the channel is considered as a 
discrete memoryless m-ary symmetric channel. These results are 
derived as functions of up/down link signal -to-noise ratio and 
the characteristics of the TDRS forward and return processors. 

5.2 Communication Link Model 

The communication system to be considered in this section is 
shown in Fig. 5.1. Information bits (assumed binary) are first 
grouped into k-bit groups and then transmitted to the satellite 
repeater by a MPSK (H=2^,k=l j2,3,, . . ) transmitter which 
transmits one of the M phase-shift-keyed' symbols per k coded bits. 
Notice that BPSK, QPSK and SQPSK signals are special cases of 
this format. The satellite repeater is a bandpass nonlinearity, 
consisting of a bandpass filter and a TWT amplifier, which in 
general has both AM to AM and AM to PM characteristics. Additive 
Gaussian noise is introduced on both the uplink (TDRS repeater 
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Fig.5:l. MPSK Satellite Communication Link 
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noise) and the downlink (ground station receiver noise). The 
demodulator is assumed to be coherent MPSK detector which 
heterodynes the received MPSK signal with a correctly phased 
local reference and decides one of M phases was sent by computing 
(in principle) the inverse tangent of the hatio between the values 
of the baseband outputs of the zonal low pass filters at the I 
and Q (inphase and quadrature) channels in the demodulator, which 
are sampled once per symbol. Fig> 5.2 illustrates the operations 
of the MPSK coherent receiver. Through an inverse mapping that 
maps each of the M-phases into a unique group of k bits the 
transmitted information can be received from the demodulator 
output phases. 

Of major concern in thissection.i*s the- effects of uplink 
and downlink noises on the performance of the MPSK communication 
system of Fig*. 5.1. Thus we have omitted the effects of inter- 
symbol interferences created by the transmit filter, the receive 
filter, and the bandpass filters in the satellite repeater. In 
other words, all the filters in Fig?. 5. 1 will be assumed zonal 
filters that pass the desired bands of signals without introducing 
intersymbol interferences, wnile limiting noise powers to finite 
values. Since for rectangular wave-shaping MPSK signals have 
equal bandwidth occupancy for equal symbol rates, this assumption 
does not hamper the tradeoff between selection of various MPSK 


systems with respect to bandwidth requirements. 


”1 


A number of papers. [5-'^-[5-^=have7'Considered the determination 
of che error probability for special cases of the channel of 
Fig: 5.1. Jain and Blachman I?-T] considered the detection of 
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BPSK signals transmitted through a hard-limited channel. 

Mizuno et-a'l.t5-2] later extended their analysis to an n-link 
channel with hard limited repeater between links. Hetrakul 'and 
Taylor [5-3], Lyons [.5-4] have also considered channels similar 
to Fig:. 5.1. There is a major difference between the approach 
undertaken in [5-T],[B-3 and that of IB-3] , [^4] . The approach of 
15-1] ,[5“3 can be called the convolution approach, which involved 
the convolution between two independent phase distributions as 
explained in detail in [5-3- The approach of [5-3] ,[5-4] first 
separated the output of the nonlinearity into a gross output 
signal term and a zero mean equivalent noise term. Through 
analytic modelling of the bandpass nonlinearity the moments 
of the equivalent noise term are then computed, from which 
the density functions of the equivalent noise terms are then 
approximated by truncated Gram-Charl ier series expansions. 

From these approximated expansions of density functions numerical 
integration are then used to obtain the probability of symbol 
errors of the MPSK system. It is felt that this approach is 
more direct. In addition to probability of symbol errors the 
transition probabilities when the channel is considered as 
M-ary channels are also evaluated. From these transition 
probabilities the probability of bit error of the uncoded 
MPSK system can be directly evaluated for any given mapping 
between the coded bits and the phase-shift-keyed symbols (for 
example, Gray codes). Also can be evaluated from these 
transition probabilities is Rq, tne computational cutoff rate, 
which is an important parameter in coding analysis. 
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A general formulation of the symbol error rate expression 
for MPSK signals transmitted through a bandpass nonlinear 
channel with both uplink and downlink noises is given in Section 5.2. 
Two different computational techniques are used to compute these 
error rate expressions, which stem from two expressions of 
the probability of symbol errors derivable from the general 
formulation. These are discussed in Sections 5.3 and 5.4. Computation 
of bit error rates and Rg from the transition probabilities are 
discussed in Section 5.5. Numerical examples for a specific nonlinear 
channel are discussed in Section 5.6. 

5.3 Probability of Error for Shuttle Signals Transmitted Through 
the TDRS Satellite Repeater 

In the uncoded case the information bits is first grouped . 
into groups of k bits and then the transmitter transmits a M-ary 
(M=2 ) PSK symbol for each k information fait. Many mappings can 
be made between the information bits and the phase-shift keyed 
symbols. One that commonly is used is the Gray code which has 
« the unique advantage that there is only one bit difference between 
adjacent phases transmitted. Let be the sequence of phases 
being transmitted, then the waveform of .the transmitter output, 
neglecting the effects of the transmit filter, can be described by 


CO 

S{t) = ^ Ap(t-kT)cos[wgt+5^j^] (5-T) 


k=-«? 


where ojg is the carrier frequency, A the signal amplitude, T 
the symbol duration, and p(t) the waveform shaping function of 
the symbols. Assuming p(t) to be rectangular pulses then the 
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transmitted signal during one symbol can be simply written as 

S(t) = A cos[(u^t+5^^] (5-2) 

for 0 _< t _< T. 

At the TDRS input, the front end noise is assumed to be white 
Gaussian with spectral density N-j . Let be the bandwidth of the 
BPF in front of ’the TWJ in the TDRS and assume B-j is wide enough 
to pass s(t) without distortion. Then the uplink signal plus 
noise input to the TWT combination' has the form 

x(t) = S{t) + n^{t) 

= (t)cos[o)^t+ni (t)] (5-3‘) 

where n-j(t) is the narrowband Gaussian noise with mean zero and 
2 

variance a-| = N-jBp and where v-j and n-j are the envelope and 
phase respectively of the narrowband process x(t). If we further 
assume the BPS's between the TWT input and at the TWT output to 
be zonal filters in the sense that they pass the fundamental zonal 
outputs of these nonlinearities without distortion while suppressing 
completely their higher zonal outputs, then the TDRS can be 
considered as a memoryless bandpass nonlinearity with AM to AM 
and AM to PM characteristics f{v^) and g(v-j) respectively. The 
output of the hDRS for input x(t) can then be written as 

y(t) = f[v^(t)]cos[o)^t + g(v^(t)) + n-[(t)] (5-4)' 

Let Ng be the spectral density of the white noise at the front 
of the receiver BPF with bandwidth 62 * and again, assume 82 to be 
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wide enough to pass y(t) wi-thout distortion, then the receive 
filter output can be represented as the following narrowband 
process 


r(t) = y(t) + H 2 (t) 

= V2(t]cos[oj^t + (5-5) 

where n 2 (t) is the narrowband Gaussian noise with zero mean and 
2 

variance 02 and where V 2 and ^2 sre envelope and phase of the 
process r(t). Since coherent detection is assumed the local 
reference of the receiver can be assumed. to track out the average 
phase shift caused by the AM/PM characteristic of the TWT so that 
the local carrier reference can be assumed to have the form 
cos[(o^t+g^, where, g* denotes the mean phase shift 

? = E[g(v-])] . (5-6) 


and where E denotes the expectation with respect to uplink-envelope 
v-|. The resulting single sample, inphase and quadrature components 
at the demodulator are then 


2 = cos[np-g] 

(5-7) 

z*s = V2 sinEng-g] 

Once every T seconds a pair of and 2 ^ is formed in the 
demodulator from which a decision is made with respect to 
which of the M phase was sent. The conventional MPSK detector, 
which is optimal for additive while noise channels, is to form 
{see Fig. 5.2) 

n2 “ g = tan‘^[z^/z^] (5-8) 
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and to decide that was the transmitted phase according to the 


decision rule: choose ^ ~ T ’ ^ ^ if and 

only if 


■'T ^ J- ir 

M " M - ^2'^ -'"mm 


(5-9.) 


for m = 0,1 ,2, . . . ,M-1 . 

Since the filtered uplink and downlink narrowband noise processes 
have uniform phase probability densities and the TORS repeater 
nonlinearity is independent of the input signal phase, there is 
complete symmetry with respect to each of the M possible transmitted 
phases. Hence without loss of generality we can consider the 
transmitted phase to be = 0. In the following we shall derive, 
in addition to the probability symbol error, the transition 
probabilities: 


VJ) =^=^1^=0 


j = 0,1,2,...,M-1 

which is the probability that the HPSK demodulator will decide 

A 

f =j /M was sent while ^^=0 was' indeed transmitted. Note that 
the probability of symbol error is simply given by 


P„(E) = l-Pj^(O) = 1-2 Pm^J^ 

J=1 


(5-ni 


which is readily available from the transition probabilities P|i^(j)- 
Recall that(eq.5-3) theinputto the satellite repeater can- be 
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written as. 


x(t) = s(t) + n-j (t) 

= V^{t)C0S[cL)^t+Tl.j(t)] 

The conditional probability density of v-j and n-j given A and ^ 
can then be readily obtained from a polar transformation of the 
joint Gaussian density function of the inphase and quadrature 


components of x(t): 


V-j 

p(v-j ,n-j jA,9^) — ^ exp 

2ira-j L 


- 2v^A cos(ri-] “'5^) 


Similarly the conditional distribution of the downlink received 


signal envelope V 2 and phase ri 2 given v-j and iT| is given by 
(see eq.5-5)-: 


2 2 

^2 r ^2 ^ ^ -2v2f(v-|)cos[ri2-n-[- 

P ( ^2 j ri2 1 ’ ni ) “ 2 2 

2,ra2 L " -202^ 


Tg(v-i)] ^ 


The conditional transition probabilities given and n-j » i '^1 >n-j ) > 

defined in an exact manner as the unconditional transition prob- 
abilities can now be obtained by integrating the conditional 

probability density p(v 2 >n 2 N-j >n-[ ) over the respective regions of h 2 
defined By eq (5-9) and the range of from 0 to 

J (2J+1)-*? 

P|v|(j iv-| ) = / I p(v 2 ,n 2 K-j >n-j )dv 2 dTi 2 (5-14) 

S(2j-1)+? *^0 

The unconditional transition probabilities Pf^(J) then obtained by 
averaging' ^1<$ over v-j and : 
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/ / P(v-),ni |A,^) f / 

■'o -'o -'o 


P(v2>ri2-N-| >n-] )dv2dri2dTi-|dv-| 
(5-15) 


In the following sections we describe two approaches that can 
be applied to further simplify the form fold integral of (5-15) in order 
to arrive at an easily computable expression of Pj^(j). 

5.4 Confluent Hypergeometric Expansion ^ Approach 
5.4.1 General Bandoass Nonl.inear Channel 

The conditional probability density of given v-j and n-j can 
be evaluated by integrating V 2 in (5-13) over the region' from 0 to «. 
Several equivalent expressions of p(n 2 1 ,n-] ) are available (see [5-8]). 
The following lists two equivalent expressions of this conditional 


probability density. 


<v2 

e ^ , r2 


P(n2lvi,ni)= cos,(n2-n2)e 


2 2- 
-p« sin (n,-no) 


f '^2 • 2 / — \ 

_ -p„ sin (no’Tip) ^ _ 

cos(Ti2-n2)® erf[p2 cos(ri2-'n2)l 


(5-16) 


p ('Ti2 1 'V-] jTi-j ) ~ 2’jj- ^ pi ^(■gi'l.) ] 2 ’^+1 ;-p2]cos n[ri2-n23 (5-17) 

n=0 


where 


^2 


^4inad pags li 


(5-18) 


^2 = n-j + g(v^ ) 


(5-19) 


In the- above equations, erf stands, for the standard error function, 
£ are the Neumann numbers 
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^ n = 0 

2 n > 1 


(5-20) 


and represents the confluent hypergeometri c function. In the 
fallowing derivations we shall use (5-1^ , which is a Fourier series 
expansion of the density of while its coefficients involve 
hypergeometric functions. 

Since 


/ 


J M M 3 


cos n[n2-iT2]d'^2 


2tT n 

-R • " ' ° 

2 mr _r. 2ir 


i — ^ + Q 

^ M M 9 


- sin cos n[j -n.|-(g(v^)~g)];n>l 

(5-21) 


the conditional transitional probabilities P^(jl ,n-j ) can be 
evaluated from (5-17) to be (note: since (5-17) is a convergent 
Fourier series, the integrals over rt^ and the summation over n 
can be interchanged in order): 


» sin 


niT 

M 


n=l 

4 iFiCf;n-H;-4]cos n[j ^ - n, -(g(v, )-?)] 
j = 0,1,2,...,M-1 . (5-22) 

For the special case when uplink noise is absent v^=A and 
0^=*? are known constants; it is clear in this’ case that the 
unconditional transition probabilities are the same as the 
conditional transition probabilities given in (5-22). In particular, 
we have, in the absence of uplink noise. 
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nir T./n.T\ 

1^2 M ^^2 ^ n p rD , 2-| 

M Tr 2-/ n r(n+l) 1^1 1^2’" 


f^(A) 


= constant 


j ^ - Cn^+g{v^)-g] 


= g(A) - g 


which is the probability of correct reception of the usual one link 

2 

MPSK system with SNR and phase incoherence li). In fact, for the 
BPSK case (M=2) it can be easily verified that 

P2(0) = ^ 

which is a well known result for coherent detection of BPSK signals. 
In general to obtain Pf^(j) we must average Pj^^(j/v-| , n-j ) over v-j 
and n-[ • Substituting (5-12),(5-l 7) into (5-15) and observing the fact 


2tt / v^A \ 2tt - 

I exp( — g- cos n-ij cos n[j - n-j -g{ v, )+g]dn-j 

\ oi / 

= cos n[j ^ - g(v,j ) + g] I^ ^ 


where are the modified Bessel functions, we have the following 

general expression for the transition probabilities for MPSK signals 
transmitted through a bandpass nonlinearity with up/down link noises: 

p (J.) . uriy 

M J TT n r(n+l ) ® 2 

0 n=l ^1 




Fi i 2^ 




j = 0,1,2,...,M-1 
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with p-j = A /2c-j , the uplink SNR, the symbol error probability is 
then given by 

P^(E) = 1 - Pj^(O) (5-27) 


Before we discuss the effects of backoff of the TWT and 
the definition of downlink SNR, we shall discuss some special 
cases of the above expressions. 

5.4,2 Reduction to the Hard Limiter Case 

In the case of the hard limited repeater, 


f(v^) 


= p 2 ~ constant, not a function of v-. 


g(v.j) = 0 (5-28) 

then (5-27) is reduced, to 

P /r-i „ M-1 2 ^ /Pi n 2-, 

' M ■ TT 2 n r(n+l) ® P2 ^’"P2^ 


n=l 




r V 

~ 


1 

"Va 


dv. 


(5-29) 


The integral in (5-29) can be evaluated by the Weber and Sonine's 
integral formula, arriving at the following result: 

_ n-iT 

M-1 2 




M 


-f z 

n=l 


sin 


M 


^rf^l 


n 


LWhTTJ 

X -jF.| E^;n+1 ;-p^] [j;n+l j-p^] 


n n 
PtP2 


picrGnsrAE p age is 

m EOOR. QUALITY 


(5-30) 


This result is in exact agreement with that obtained by Mizuno et al . 
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f P(v,) 



where p(v-j), defined as 


7?b: 




X cos n[g(v-|)-T|)])dv-j 


P(v-j) 


“2 


is readily recognized as the probability density of the uplink signal 

2 

plus noise envelope. Clearly when the uplink SNR p^-x=» (i.e. when 
a-|->0) we have 

lim p(v,) = 5{vn-A) (5-3 

where 5 is the Dirac-Delta function; and also we have 


lim — > - n — V - = 1 

■’*" '"(it) 

Hence the integral (5-34) becomes, as cr-j-»-0, the following 


Pg -|F-][^;n+l ;-ppcos n[g(A)-g] 


Substitute (5-39) into (5-27), we have 


11m PJE) = 


a.j->0 


2 ^ will 

IT 2^ n 
n=l 


r(-^l ) 


r(n+l ) ^2 1*^1 ^2’ ’“^2^ (5-40) 


X cos n[g(A)-g] 
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where ^ f(A)/v^a 2 in this case defines the dov/nlink SNR. Here (5-40) 
is [see (5-23)] clearly the probability of error of a one- 
link MPSK system with phase incoherence g(A)-g’, as expected. 

Next we consider the case when the downlink SNR 
(i.e., a2->0). Using the known relationship on confluent 
hy pergeometri c functions 

^F^[a;3;x] = e^ ^F^[e-a;a;-x] (5-41) 

and the asymptotic behavior 

eV®; x » ] (5-42) 


eq. (5-27) can be reduced to 



/Av. \ 

X ^ny~2 (5-43) 

With the following relationship, which can be easily verified by 
using the Jacobi-Anger formula; 




j 2i ) cos “ ff ^ n 


/A'^-j \ _ 

X ^jcos n[g(v.|)-g] 


(-5-44) 


the asymptotic P given in (5-43) is verified to be in agreement 
s 


with the following. 


ORIGINAU PA:GE5 IS 
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lim P„(E) 



f -g(v^)+g 

J -g(v-,)+? 


Jl. 

2irJ^ 


exp - 


'A +v-j“2Av-| cos n-j 


2a 


dn-jdv| 


(5-45) 


which is easily seen to be a correct result when downlink noise 
is absent, from basic arguments. 

5.5 Backoff and Downlink SNR 

Unlike the case of the hard limiter channel, the probability 
of symbol error expression ^-27) for a. general bandpass nonlinearity 
is not explicitly .given in terms of up/down link SNR's. In fact, 
both and a^ are present in (5-27) . It is useful to normalize (5-26), 
(5-27) in terms of SNR's and certain parameters of the 
nonlinearity which we will define in this section. First 
let f^ denote the AM/AM characteristic of the limiter and let 
v-j denote the input envelope at which the limiter saturates. 
Correspondingly, define the limiter backoff coefficient to be 
the rati o : 


(5-46) 


"1 ,SAT 


2 

Then, in terms of and v.j uplink noise variance a.j in 

(5-26), (5-27) can be written as; 


2p 


^ L^l ,SAT 

, 2 
2pT 


(5-47) 


Secondly, let f-j-jg-p respectively denote the AM/AM, AM/PM characteristics 
of the TWT, and P.^ denotes the input power at which the TWT 

saturates. The backoff of the TWT, By, is customarily referred to 
as the ratio between the input power into the TWT and P^^ 

P. 
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in, SAT 

This backoff is controlled by the adjustable drive between the 
limiter and the TWT as illustrated in Fig. 5. 1. Let a be the 
gain of this adjustable drive. Then a is related to through 
the relationship: 


a = 


pSAT 

^in 


r^)/z 



where fj^(v^}/2 is the average power of the limiter output: 


(5-48) 


T _ 1 r -2, V ^1 Z ^ 

2 “ 2 I \^^1^ 2 ■ jPi 

n V 


2 , '’1 


1 2 a?rn 


IJ —k dv. 


In terms of py and v.j eq. (5-50) can be written* as 

( p~\ \ 'i'2 

® L '’l ,SAT \ r 2,. . ^1 ,1^ u ^ ^ 

j~Z ^1 ) ~2 


where we have defined the normalized envelope: 

un i/Tp, , 


— ^1 


^ L ^1 ,SAT 


The effective AM/AM and AM/ PM characteristics f and g of the 
combined nonlinearity of the soft limiter and the TWT with 
their respective backoffs are then given by: 


9(^,) = y,); 


' rJLinC^i 
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where a is given by eqs. (5-49) and (5-51). 

Another parameter that needs to be defined in the error rate 
performance is a measure of the downlink SNR. It is felt that in 
the case of a general bandpass nonlinearity it is most straight- 
forward to define the downlink SNR at saturation (of the TWT): 


2 _ ^0UT,SAT (5-54) 

^2, SAT I 


where is the output power of the TWT at saturation. Notice 

this definition is consistent with that of the downlink SNR in 

2 

the hard limiter case. With this definition of sAT 

together with 3j^(5*46), B.j.(5-48),a (5-49), f and g^-53), the transition 

probabilities Pjv^{j) (5-26) can be written as 





r(^l) 
n r(n+iy 


2 

v-| exp- p-| 



X Ij^(i/rpi v-j ) 


' P 

7 ^ 


2, SAT 
OUT, SAT 






2, SAT 
OUT, SAT ^ 



2 

X cos n[j - g(v 


) + g] > dv. 


1 


j = 0,1,2,... ,M-1 (5-55) 


in ^-55) the functions f and g are those defined in (5-53), which 
are dependent upon the backoff coefficients and Bj of the 
limiter and the TVIT respectively, in addition to their respective 
AM/AM and AM/PM characteristics. Eqs. (5-55) gives an expression of 
the error rate which is easier to apply when performing tradeoff 
between various up/down link SNR's and backoff conditions. 


-MQ- 


InL^m 



■J!tnCc 


om 


5.6 Error Function Expansion Approach 

An alternate expression of the symbol error probabilities can 
be obtained by first writing the conditional probability of symbol 
errors given uplink envelope {‘v.j) and phase (n^) in terms of 
error functions. Recall from eqs, (5-14) and (5-1^ , the probabilities 
of correct detection, Pj^(O), are obtained by averaging the 
conditional correct detection probabilities (see eq. (5-14)) 


IT f „ 

M + g 


P{^(0|v-|>ni) ~ 


■ / /; 


^2 * ^2 ^ ^1 * ^1 (5—56 ) 
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where the conditional density p(v,;i 2 |v, .n-, ) is given in (5-13) 
Defining the new variables x = and letting f(v-j) = f(v.j)/a 2 , 

then (5-56) can be written as 
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x^+?^(v.j )-2xf(v.| )cos[Ti 2 -ni-g(v.j )]■ 

(5-57) 
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X cos[n 2 "n-|-g(v-j )] = x cos TI 2 cos[iTj+g(v.j )] + x sin ri 2 sin[n-j+g(v.j )] 

(5-58) 

and introducing the transformations z = x cos n-j and u = x sin n-| 
into (5-57) we can write 
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z tan ^ 


-z tan ^ 
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(z-V^+(u+0)^ '' 
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dto 

(5-59) 


where the I-Q channel components are defined as 
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I = f(v^ )cos[Ti^+g(v^ )+g] 

^ — 

q = -f(v-])sin[n-|-f-g(v-j)+gl 

With X-j = Z-I. and ^2' = u+Q eq. (5-57) can be simplified to 

Pn^(Qlv^ ,n-]) = f exp(-x^/ 2 ) C exp(-A2/2)dX2 

•'-I •'a 


where 


a = -[(a-j+I) tan ^ - Q] 
b = C(A-[+I)tan Q] 


Since 


erf(x) = -erf(-ix) 


X .2 


erf(x) = 1 - erfc(x) - J’ e"^ 


we can rewrite (5-61) as 
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/(A,+I)tan J -q \\ 

+ erf f — ! ) 1 dA. 


This is a general result which we proceed to simplify. Let 
z = x+I then (5-58) can be written as 

\ 

PmCoNv^t) =• ^ J f(i.Q,x)dx + f(i,q,-i) 


where 
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(x+I)tan ^ +Q 

f{I,Q,x) ^ expC-x^/2) / exp(-yV2}dy 

-[(x+I)tan ^ -Q] 


Now f(I,Q,-I) = 0 so that 
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To evaluate this integral we need to evaluate sf/sl; thus 
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which can be simplified to 


It = tan i^exp 
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where 


a ^ I sin^ g + Q sin | cos ^ 


b ^ I sin J q cos J 


M 


c ^ I sin^ ^ - Q cos ^ cos ^ 


d - I sin ^ - Q cos ^ 


M 


M 


This simplifies 3f/3l as much as we need for the present. 

We now continue with thedevelopment. Substitution of (5-69) 


(5-66) 


(5-67) 


~ -01 
(5-68) 


(5-69: 


(5-70 


into (5-67) and making the change of variables 



leads to 




(5-72) 


after considerable manipulation and simplification. The conditional 
error probability is then found from the above expression by 
integrating between 0 and I, i.e., 

I _ (5-73) 


Since 1=0 implies f(v^) = 0, then 




(5-74) 


and we have 




I 9Pf^(0lv.^ ,n-j) 


The conditional probability of symbol error given and n-j is 
then given by 


P„(E|v,,ni) = 


-I 


9Pjv|{0|v.| ,n-|) 


when (5-72) is substituted into the integral of (5^75), we obtain 
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■f ( V, ) _ _ 

^ 7 2 ^~ Sin[ ^ +(n^ + g(v^) + g)] 

The overall probability of symbol error can then be evaluated 
by averaging ^-77} over v-j and , as indicated in eq. (5-15). Some 
simplification of (5-77) can be made for the cases of M=2-and M=4. 
Their results can be shown to be the following. 

BPSK (M=2) 

, f(v,)cos[n, + g{v]) + 

P2(Elv^,n^) = l-erfc 


QPSK (M=4) 
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P^(Elv.j,n-}) = ^ [erfc(A^) + erfc(A_)] 


- 1 erfc^(A^) - ^ erfc^(A^) 


The same approach as that of Section 5.3 can be used here to 
take into account the limiter and TWT backoffs and up/down link 


a 






PROBABILITY OF ,ERRO£. 



DOWN-LINK SNR (DB) 


Fig^.3. Symbol Error Probability vs Downlink SNR with Uplink SNR as 
a Parameter. 
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Fig. 5.4 Symbol Error Probability vs Downlink SNR with Uplink SNR 
as a Parameter. 
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6.0 ROCKET EXHAUST EFFECTS ON SHUTTLE RF COMMUNICATIONS 


In Vw- <om 


6.1 Introduction 

The Shuttle Space Transportation System' has two Solid Rocket 
Boosters (SRB) attached to the main propulsion unit. The two SRB's 
assist the main propulsion unit during liftoff and into the early 
and critical phase of the launch trajectory, depending on the 
look angle from the ground station(s) to the Shuttle antenna, since 
SRB plumes are relatively opaque to RF signals. The uncertainty of 
the efforts of the SRB plumes on Shuttle communications is further 
compounded by the fact that these SRB's are the largest ever used. 

Thus, there is no directly applicable past experience that can be 
depended on explicitly. 

Certain extenuating factors make it difficult to predict the 
plume effects on the Shuttle RF communications. These factors 
include ground reflections of the RF signals, the cloud of ground 
debris stirred up by the wind and pressure generated by the rocket, 
dynamic effects of changing altitude and look angles. 

Knowledge of the effects of the TRIDENT I rocket exhaust on RF 

transmission are perhaps helpful to run some RF signal tests for predic- 

ting Shuttle exhuast (plume) effects on the RF transmission. The 
RF tests can be designed to minimize these extenuating factors and 
maximize the TRIDENT test results can be extrapolated or help in 
designing tests for determining SRB plume effects on the Shuttle 
Communications during launch. 

The Shuttle Communication frequencies at S-band are: 

*2287.5 MHz PM downlink TLM 

*2106.4 MHz PM downlink CMD 

•2217.5 MHz PM downlink 
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This report describes the measured results and the effects of 

the TRIDENT I rocket exhaust (plume) on the amplitude and the 

phase of the radio frequency transmission at the GPS/SATRACK 
frequency of 1575 MHz. The contents of this report and the 

measurement results are based on the tests conducted by Lockheed 
Missiles and Space Company and documented in the report entitled, 
"TRIDENT I Rocket Exhaust RF Transmission Measurement Report."* 

The measurements were made before and during the following motor 
firing at the Naval Weapons Center, China Lake, California; 


CST 1 Third Stage 
CST 2 Second Stage 
CST 3 First Stage 


The following conclusions were reached. 

1. The three firing tests successfully measured attenuation, 
phase shift, attenuation jitter, and phase jitter caused 
by the rocket motor exhaust. 

2. Empirical formulas quantifying the relation between attenuation 
and both jitters are' developed. The relation can be expressed 
as follows: 

Attenuation Jitter(dB) ~ 0.4xAttanuation (dB) 
Phase Jitter (degrees) ^ A.OxAttenuation (dB) 

3. The simulation of the rocket motor exhaust by aluminum foil 
strips prior to each firing aided analysis of the results. 

4. The correlation between field measurements and the model 
are considered good enough to warrant the use of the model 
in planning future tests. 


*The contents of this memo are extracted from the following report. 
"TRIDENT I Rocket Exhaust RF Transmission Measurement Report," 
Lockheed Missile and Space Company, Inc.; Report No. LMSC-D057671 , 
March 1977. 
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Section 6.2 of this report contains a description of the test 
equipment and the test installation. 

Section 6.3 of this report contains the observations made 
during the motor firings, during the simulation of the plume 
by aluminium foil strips, and by runs of the model. 

Section 6.4 of this report analyzes the correlation of the 
observations within each set and between sets while the appendix 
presents recorded data. 

6.2 Test Description 

Three combined systems test (CST) test firings were conducted 
and two RF transmission systems were used. Figure 6. 1 shows a 
functional block diagram of RF transmission system for test Setup 
No. 1. An elevation sketch of the antenna configuration is 
illustrated in Figure 6. 2. During CST No. 1, the RF transmission 
system consisted of a CW oscillator feeding a fixed antenna located on 
one side of the rocket plume,, with a vertically movable receiving 
antenna located on the opposite side of the plume. The oscillator 
frequency was 1575 MHz. A network analyzer was used to compare a 
reference sample of this transmitted RF signal to the signal 
received by the moving antenna. The resulting phase and amplitude 
output signals were recorded on a magnetic tape. 

Figure 6.3 shows test configuration of CST Tests Nos. 2 and 3. 

For -the CST tests 2 and 3, a modified system was employed which 
provided the capability of transmitting from an alternate antenna 
placed approximately 150 feet from the plume or from the 
original transmitting antenna which was located approximately 
25 feet from the plume. Also, when transmitting from the aft 
antenna, the forward antenna was connected to an auxiliary 
receiver in order to measure the RF power reflected from 
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Figure 6.1. Rocket Motor Plume RF Transmission Test Set-up No. 1 
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Figure 6.2. Rocket Plume RF Transmission Test Evaluation Sketch. 
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Figure 6,3. Rocket Plume RF Transmission Test Set-up No. 2. 
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the plume or from the aluminium foil plume simulation surface’. 

The simulation experiment was designed to approximate plumes of 
various widths. 

6.3 Test Results 

6.3.1 Firing Results - Forward Antenna 

The observations from the three firing runs were reduced and 
plotted in Figures 6.4 through 6.6. The four observations on one 
plot permit the relations of the four elements (phase, attenuation, 
phase jitter and attenuation jitter) to be compared. 

6.3.2 Simulation Results - Forward Antenna 

Figures 6.7 through 6. T2'. summarize the observations. The 
simulation experiments were performed before each firing and are named 
for the firings they preceded. 

6.3.3 Model Results - Forward Antenna 

A model to forecast the effect of rocket plume on the RF 
transmission has been detailed in a Lockheed Missiles Space-Company 
report entitled, "C4 Rocket Exhaust RF Transmission Model," LMSC/ 
D057575, 8 December 1976. Figures 6.13 and 6.14 show model 
predictions for the forward antenna configuration. 

6.3.4 Firing Results - Aft Antenna 

Figure 6.15 shows the observations of the first stage firing 
from the aft antenna. 

6.3.5 Simulation. Results - Aft Antenna 

Figures 6.16 through 6.19 show the results of the two aft 
antenna simulation observations. 

6.3.6 Model Results - Aft Antenna 

Figures 6.20 and 6,21 show the results of runs of the model with 
different plume widths^ for the aft antenna configuration. 
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Figure 6.4. First Stage Firing Summanv, Forward Antenna. 
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Figure 6.5 Second Stage Firing Summary > Forward Antenna, 
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Figure 6.6 Thrid Stage Firing Summary, Forward Antenna. 
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Figure 6.7 First Stage Simulation, Attenuation, Forward Antenna. 







Figure 6.-9 Third Stage Simulation, Attenuation , Forward Antenna. 
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Figure 6.11 Second Stage Simulation, Phase Shift, Forward Antenna. 
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Figuy'e6.l2 Third Stage Simulation, Phase Shift, Forward Antenna. 
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Figure 6,13 Model Runs, Attenuation, Forward Antenna. 
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Figure 6.14 Model Runs, Phase Shift, Forward Antenna. 
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Figure 6.16 First Stage Simulation, Attenuation, Aft Antenna. 
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Figure 6.18 First Stage Simulation, Phase Shift, Aft Antenna. 
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Figure 5-19 Second Stage Simulation, Phase Shift, Aft Antenna. 
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Figure 6.20 Model Runs, Attenuation, Aft Antenna. 







6.4 Analysis - Correlation of the Observations 
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6.4.1 Firing Correlation 

Figs. 6.22 and ’6.23 show plots of the forward antenna attenuation 
and phase shift of the three firings. The attenuation comparison 
shows all three stages converging to zero at about 10 feet (antenna 
height) but they show different magnitudes and shapes near their 
maximum values. Similarly, the firing phase comparison shows fair 
correlation for the maximum values. 

6.4.2 Firing Jitters 

Figures 6.4 through 6.6 show that the attenuation and phase 
jitters have a rough correlation to the magnitude of the attenuation 
for the largest values of attenuation. Note that the jitters do 
not vanish when attenuation goes through zero. The amount of 
mechanical jitter present in the data and its contribution is 
unknown. 

Based on the test results and correlation of the observations, 
the peak values of attenuation, attenuation jitter, and phase jitter 
are tabulated below in Table 6. 1. 

Table 6.1. Peak Values 



Parameters 

Stages 

Sum 

IFWD 

TAFT 


3 

1 

Attenuation (dB) 

20 

17 

20 

15 

72 

2 

Attenuation Jitter 

6 

5 

10 

9 

30 


(dB) 






3 

Phase Jitter (dB) 

40 

32 

120 

80 

272 


Attenuation Jitter/Attenuation = 30/72 ~0-4 dB/dB 
Phase Jitter/ Attenuation = 272/72 /n- 4 degrees/dB 
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Figure 6.22 Firing Comparison, Attenuation, Forward Antenna. 





From the above it can be summarized that the following emperical 
formulas quantifying the relation between attenuation and both 
jitters can be used. Note that the jitter values are peak-to-peak 

Attenuation Jitter in dB 0.4 (Attenuation in dB) 

Phase Jitter in degrees ^ 4 (Attenuation in dB) 

6. 4.3 Simulation Correlation 

The simulations repeated the same experiment and gave similar 
results within the limits of experimental error. Figures 6.24 and 6.25 
show the case of four feet width foil in the forward antenna set. In 
the forward antenna set the two and four feet width foil results 
had a good correlation. The correlation between the runs decreased 
as the width increased. The aft antenna set had similar results 
as forward antenna case, i.e., as the width of foil increased, the 
correlation between the runs decreased. 

6.4.4 Triple Comparison 

The forward antenna data (attenuation and phase comparisons) 
are plotted in Figs. 6,26 and 6.27 for a four feet comparison, 
to correlate and see how closely the three kinds of data compare 
(i.e., firing third stage, simulation and model). Similarly, the 
aft antenna observations are plotted in Figures 6,28 and 6.29 as 
a six feet comparison. 

The simulation data show the most rapid swings, the firing 
data changes less abruptly, and the model results are smoother 
and shallower than either. The simulations provided a spectrum 
of widths to roughly approximate the plume. 


ORIQINAB EAGE IS 
OB EOOK 




IGHT (FEET) 


00 

00 

1 



ATTENUATION (DD) 

Figure6.24 Simulation Comparison, Attenuation, Forward Antenna. 
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Figure 6.25 Simulation Comparison, Phase Shift, Forward Antenna, 
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Figure 6i26 Attenuation Comparison, Forward Antenna. 



r 



Figure 6.27 Phase 
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Figure 6.29 Phase Shift Comparison, Aft Antenna. 
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RECORDED DATA 


F1g, 6 A.-l Illustrates the analog recording of the firing 
data of the second stage using twenty-five cycle filters on the 
amplitude and phase data; runs using this filter bandwidth 
were used for the basic results. Fig. 6 A. 2 represents a 
record of the second stage firing data using 5 kHz filters. 

The peak-to-peak deviations on runs with this filter were 
measured and reported as phase and amplitude jitter. 

As the antenna is cycled up and down with no obstruction 
between transmitter and receiver, the received amplitude and 
phase both change. The amplitude changes are caused primarily by 
the interaction of the beams of the two directional antennas used. 

The phase changes reflect the fact that the distance between the 
antennas is changing. To remove these changes caused by variable 
antenna position, phase and amplitude measurements made during a refer- 
ence run with no obstruction are subtracted from those with an 
obstruction and the difference attributed to the obstruction. The 
primary indication that data is valid is the fact that it approaches 
free space, i.e., converges to zero difference, above the obstruction 
whose centerline is at zero feet for all tests. 
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Figure 6A.1 Analog Record 



AMPLITUDE 




fi .,kf 







\ 




li !i I I 






ANTENNA HEIGHT 


44, t t tit t 4 4tt 4 • t 4 *4 t i tit 4 4 4«t 4 » 4 4 4 444 4 4 444 4 «t i 

“ * i i «i * t i tut 4«tl t Mi t < it t ft ittt t«*t tft »t f 'i t t >f » » t Mtt tttt 

|l*tf!rtl»f*f'*.»f tf**ff»i«»W|f«fltff t m titn nm yBiigi»HHiHt»tM»^fHatMM{HHt4*i»4ttWfWf<H»»« 

Figure ^A.‘2 Analog Record. 




7.0 PERFORMANCE SENSITIVITY ANALYSIS TO DATA ASYMMETRY FOR 
UNCODED AND CODED SHUTTLE BPSK KU-BAND RETURN LINK SIGNALS 

7.1 Introduction 

The performance sensitivity to data asymmetry is analyzed for 
NRZ data, both coded and uncoded. The definition of data asymmetry 
in agreement with the TDRSS Users' Guide (STDN No. 101.2, Revisions) 
The math model used and assumptions made are summarized. The CNR 

degradation is tabulated for the infinite bandwidth linear white 

-2 -3 

Gaussian noise channel for bit error probabilities of 10 ,10 , 

10 lO"^, 10~^ and data asymmetries of 5, 10, 15, 20, and 25%. 

In the uncoded-BPSK case, both periodic and random data is 
considered. Degradations for the K=7 rate one-half convolutional 
code is tabulated over the above range of data asymmetry and 
bit error probabilities. A more elaborate analysis of the coded 
case is documented which takes the effect of filtering into account 

7.2 Data Asymmetry Degradation on the Linear White Gaussian 
Noise Channel 

7.2.1 Definition of Data Asymmetry 


For 100% transition density the data asymmetry can be defined 
as (Ref. 7-1) 


"f" •• 

longer bit - shorter bit ^ T -T 
longer bit + shorter bit T^+T' 

— + 

I 



Figure 7.1 Illustrates Asymmetric Pulse Stream. 
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X IS given in Table 7.1 for the 3 three-bit combinations which 
are affected by the data asymmetry. In the remaining 5 cases 
it is 

X = AT ■ (7-4) 
Note that in Table 7.1 optimum timing was assumed for the detector. 
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From (7t2), (7-3), (7-4) and Table 7.1 the average error 
probability follows easily; 

1. NRZ random bit stream (50% transition density) 



2, 100% transition density' (periodic data) 
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Table 7.2a. Degradation in dB for Uncoded Signals (Periodic Data). 


El 

5 % 

10% 

15% 

20% 

25% 

10'^ 

.5 

1.2 

2.1 

3.4 

5.0 

10'^ 

.5 

1.4 

2.5 

3.8 

5.4 

10-'^ 

.6 

1.5 

2.7 


5.6 

10“^ 

.6 

1.6 

2.8 

El 

5.7 

10^^ 

.7 

1.7 

2.9 


5.8 


Table 7.2b. Degradation in dB for Uncoded Signals (Random Data). 


Bi 

5% 


15% 

20% 


10"^ 

.2 

.6 

1.1' 

1 .7 

2.7 ' 

10'^ 

.2 

.7 

'1.4 

2.5 

3.9 

lo-'t 

.3 

.9 

1.8 

3.0 

4.6 

10'^ 

.3 

1.0 

2.1 

3.4 

5.0 

10"^ 

.4 

1.2 

2.3 

3.6 

5.2 
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2 

where Ej^ = A T. The necessary increase in E|^/Nq to offset the error 
rate degradation is tabulated in Table 7.2a for periodic data and 
Table 7.2b for random data. 

7.2.4 Error Probability for Coded Signal 

We consider two models, one which includes the effects of 
filtering and leads to bit-error bounds and one based on an infinite 
bandwidth channel. For the latter, performance results are tabulated. 
7.3 Description of Bandlimited System 

A communication system model is shov/n below. The memoryless binary 
source produces a stream of binary bits, each equal to 0 or 1 with 
equal probability. The convolutional code chosen is one that will 
be used on the Space Shuttle and is K=7, rate = and nonsystematic. 
The modulator creates a baseband BPSK signal corresponding to the 
coded bit stream, a sequence of pulses with amplitude either A or 
-A for some A. Asymmetry of- the pulses is allowed; specifically. 



Channel 


Detector 
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transitions from A to -A occur late. The channel consists of a 
Chebyshev filter specified by its degree n and bandwidth B, followed 
by the addition of white Gaussian noise. A signal -to-noise ratio 
is specified. The Viterbi decoder used here ignores intersymbol 
interference. It truncates beyond 4K branches back. 

7.4 Justification for the Detector 

In order to partially justify the detector we have chosen, 
we first considei^ the simpler system* 



Let us say that the binary source emits a total string of length N; 
then the length of the coded string is 2N. Suppose Cj^ = +1 , 
k = Q,...,2N“1. The modulator produces the signal 


s(t) 


where 




pjt) 


0 ^. 


2N-1 

Z ^kPk(^) 

k=0 

A for k j £ t < ( k+1 ) |- 
0 otherwise 
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The set of functions pj^{t)/k = 0,,..,2N-1, defined on 0 ^ t ^ NT 
is an orthogonal set of equal -energy signals, when the inner 
product (x,y) of two functions x(t) and y(t) is given by 

•NT 


(x,y) = 


= / 

J n 


x{t)y(t)dt 


Let r(t) = s(t) + n{t) where n(t) is white Gaussian noise, 

•1 

Let {s (t): i} be the set of all possible modulated coded 

bit strings. Then we define the vectors _r with components 

i . , 

s^, given by 


sj, = (s\pj^) 

k = 0,...,2N-1, all i 

"k = ^^’Pk^ • 


Then [1 ,§2.1 ,§2.2] the maximum-likelihood decision rule is to 
decide that s^(t) was the transmitted message when i maximizes 


2N-1 

Z 


k=0 



This, is what is done in the detector above. The integrate- 
and-dump part produces the vector _r, within multiplication 
by a constant; the Viterbi decoder takes the vector inner 
products of r.with all possible s_^ and yields the data string 
corresponding to the s^ that gave maximum inner product.' 

The decoder in our complete system is not optimal. The 
filter in the channel produces intersymbol interference, which 
our Viterbi decoder will ignore. Also, bit asymmetry throws 
off decoder performance. A Viterbi decoder can be devised which 
accounts for intersymbol interference [3, §4. 10]. 
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7.5 Bit Error Bound 

The following development of a bit error bound is based on 
results given in [4]. The modulated signal s(t) will consist of 
a string of 2N rectangular pulses, each of duration T/2. We allow 

bit asymmetry, where 
asymmetry = n 

_ duration of longest pulse-duration of shortest pulse 
duration of longest pul se+duration of shortest pulse 


This can be modeled by having transitions from 1 to -1 occur 

T 

late by a time of A sample piece of signal is shown below. 


1 -■ 
0 - 



We define two "pulses" p-|(t), P 2 (t), shown below. Now we can 
express s(t) as a sum of pulses. Let P be the baseband signal 


P^(t) 


T 

2 


T I — 

02 j 


P2(t) 
t 


power, so that the energy per data bit is PT. Let Cj^, k=0, . . . ,2N-1 , 
be the sequence of code bits, Cj^ = +1 . Then 


2N-1 


s(t) = /r 2 j ) 

k=0 


OM51NAI3 ESGSEI 13] 
m KK)K QlSEEESj . 
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where 


Pk-i.kf^^ 


P2(t) if = 1 and Cj^ = -1 
, p-j(t) otherwise 


and c ^ s -1 . 

We now pass s(t) through the filter, which has impulse 
response h(t). Let 

s' = s«h 


pi = p.ah , i = 1,2 

Pj_] 1 Pi_i ^ " 0,...,2N-1 . 

Then 

2N-1 -r 

s'(t) = Y. '=kPk-i r > 

k-0 


To s*(t) we add noise n(t), which has two-sided p ower 
spectral density Nq/ 2, obtaining 

r(t) = s' (t) + n(t) . 


Now the signal enters the demodulator. First the signal is 
delayed by -tg in order to compensate for the delay caused by the 
filter. Then the vector _r is obtained, with components r^, given by 

r{k+l)T/2 

^ 'VT J r{t+tQ)dt, k = 0,...,2N-l 

kT/2 


The vector r is the sum of two vectors s' and n, where 



jCinCL 


om 


rp AM]J/Z 

k " #J 

kT/2 

Y r(k+l)T/2 
kJ/Z 




s ' (t+tQ)dt 


n(t+tQ)dt 


Consider the probability of an error event. The vector £ 
is the coded bit string that was sent. Let £ be another one of 
Hamming distance d from £ then the pairwise error probability that 
£ is chosen over £ is 

= Pr[(r,c) > (r,c)] 

= Pr[0 > (£,£-£)] 

/V A 

= Pr [-(£,£-£) > (£',£-£)] 


First, we look at the term -(£,£-£). It is a zero-mean 
Gaussian random variable with variance given by 

2M-1 


A 0 ^ilZJ 9 *^9 

EC(n.£-c)^] = E y" ^(c.-cj2 


k=0 

2N-1 , . , 

E 

k=0 


We find that 


En 


f E 


■(k+l)T/2 


n(t+tQ)dt 


'(k+l)T/2 


n(T+tQ)dx 


'kT/2 

2 


'kT/2 


dt 


'kT/2 


N 
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Nq 

So the variance of “(n.*£~£) is d* "45 or 2dNQ. Therefore 


om 


' >c-c ) ^ 


where 




Q(x) = /^ I 


t)h 

Now we look more closely at s'. Its k~ component is 





We define, for i = 1 ,2 


p'(t+tQ)dt, k= [- ,...,2N-1 



We assume 9^^^ ~ 0 for k < -u and k > v, t = 1,2. The delay -tg 


was chosen so that is the largest of all i=l ,2. 


We 


define 


^ 1 - 1,1 




Si 


(k-i) 


If = 1. c, = -1 


otherwise 


Then 
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om 


■i - ^ Z 

i=0 


k+i 


= ^ H s-4-1.1 

i=k-v 


Now we can bound (^‘ ,£-£). 


2N-1 

(5'.c-c) = J2 
k=0 

2N-1 k+. 


E E <^i4^\i(v^k> 

k =0 i=k-v 
^ N -1 / . ^ k+y f ,\ 

- E'^tVk-i.k" E '=i3'-i.i>'v=k) 

k=0 i=k-v 

- /F V c.'^) Vc.<^) (c a) 

Xu k^k-l ,k^^k ^k^ Lu Z_. ^i^i-Ui^'^k’^k^ 


k=0 


k=0 itk 


^ /P 92°^2d - yF2d ^ max(jg|^*^[,|g^^'^|) 

1=-y 

i70 

= /F2d[g^°^ - ^ max(|g|'^[,[g^'^|)] 

i=-y 

i?^0 


Therefore 




g<“> - max(|g'’>|,|g<')|) 


1— y 

lYO 


The values of oj* jg^ can be found for any given filter and amount 
of asymmetry. For example, for a fifth-degree Chebychev filter 
with bandwidth equal to four times the data bit rate and ripple of 

OEIQINAI5 PAGE IS 
OE P-OOP, QUAPyCS 
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V 

10%, asymmetry of 0 gave | , |gl^'^ | ) 

i=-U 


ijiO 

a value of .58, while asymmetry of 25% gave a value of .32. 
The bit error probability is bounded by 


JZinCZ 


1 

om 


’’b ^ E E 

n-l 

where a(n,d) is the number of paths of coded weight d with n data 
bits different from the transmitted data bits [2, §4. 4]. For high 
signal-to-noise ratios this is dominated by the minimum-distance 
path, so that 


P 


b 


< 



e 


where e is small compared to P. . For the code chosen, d „ = 10. 

7.6 Performance Degradation for InfTnfte Bandwidth Channel Case 

Using the transfer function bound the performance can be computed 
for a Viterbi decoder. The channel can then be represented-by a 
model for the transition probabilities given by 



Using the BPSK error probability 



= j (p+q) 
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ont 


the performance for the coded signal can be found using standard 
conversions of to The CNR degradations for a coded channel 
are given in Table 7.3 for a rate one-half K=7 convolutional code. 


Table 7.3. Degradation in dB for Coded Signals. 
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Figure 7.2. Channel Model. 
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